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Abstract 



Tremendous progress has been achieved in neutrino oscillation physics during the 
last few years. However, the smallness of the #13 neutrino mixing angle still remains 
enigmatic. The current best constraint comes from the CHOOZ reactor neutrino 
experiment sin 2 (26 13) < 0.2 (at 90% C.L., for Am 2 tm = 2.010" 3 eV 2 ). We pro- 
pose a new experiment on the same site, Double-CHOOZ, to explore the range of 
sin 2 (2613) horn 0.2 to 0.03, within three years of data taking. The improvement 
of the CHOOZ result requires an increase in the statistics, a reduction of the sys- 
tematic error below one percent, and a careful control of the cosmic ray induced 
background. Therefore, Double-CHOOZ will use two identical detectors, one at 
~150 m and another at 1.05 km distance from the nuclear cores. The plan is to 
start data taking with two detectors in 2008, and to reach a sensitivity for sin 2 (2613) 
of 0.05 in 2009, and 0.03 in 2011. 
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Chapter 1 

Physics opportunity 



Neutrinos play a crucial role in fundamental particle physics and have a huge 
impact in astroparticle physics and cosmology. Before 2002, neutrino oscillation 
physics was still in a discovery phase, even though strong evidence for atmospheric 
|iMB91l ISUTT971 ISKMl ITyT&HMl IRonOll ISK02aj and solar neutrino oscillations have 
already been established since 1998. Thirty years after the discovery of the so- 
lar neutrino anomaly |Cle98l ISAG02I IGAL99MGNO00|, the com bined SNO Super- 
Kamiokande discovery of the flavor conversion SK0213 ISNU02| together with the 
first reactor 77 e flux suppression observed by KamLAND KAM02 , is now moving 
neutrino physics to a new era of precision measurements. 

In the Standard Model of electroweak interactions, neutrinos are massless par- 
ticles, and there is no mixing between the leptons. There exists only a left-handed 
neutrino, and a right-handed antineutrino. In the quark sector of the Model, the 
mixing between quark weak and mass eigenstates occurs among the three flavor 
families, and the amount of mixing is determined by the CKM mixing matrix. In 
the lepton sector, the analogue of the CKM matrix for quarks is just the identity 
matrix, and three conservation laws have been empirically included, for the three 
lepton families. 

The strong evidence for non-zero neutrino masses clearly indicates the existence 
of physics beyond the minimal Standard Model. The smallncss of neutrino masses 
together with the amounts of lepton flavor violation found in neutrino oscillation 
experiments provide insights into possible modifications of the current Standard 
Model of electroweak interactions, and open a new window towards the Grand Uni- 
fication energy scale |Cel79| . 

In the current paradigm, the neutrino mass and weak eigenstates are related 
through the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) neutrino mixing matrix 
|Mak621 IPon58| . A synthesis of atmospheric, solar, and reactor neutrino oscilla- 
tion data requires the existence of (at least) three-neutrino mixing. The PMNS 
mixing matrix can be parameterized by three mixing angles so i, 6* a t m , #13, and 
one or three CP-violating phases, depending on the Dirac or Majorana nature of 
the massive neutrinos |PDG00j . Although not favored by the current data, a sce- 
nario with more than three neutrinos might be required to account for the LSND 
anomaly LSN98 . In this case, the mixing of the three active neutrinos with the ad- 
ditional sterile neutrino(s) decouples from the oscillations described by the PMNS 
matrix. The presently running MiniBoonc experiment will settle the controversy in 
the near future |BOO02| . 



11 



12 



CHAPTER 1. PHYSICS OPPORTUNITY 



atmospheric + K2K : 




Figure 1.1: Solar and atmospheric allowed regions from the global oscillation data 
analysis at 90 %, 95 %, 99 %, and 3<j C.L. for 2 degrees of freedom |Mal03| . 



A wide range of experiments using accelerator, atmospheric, reactor, and solar 
neutrinos will be necessary to achieve a full understanding of the neutrino mixing 
matrix. 

Solar neutrino experiments combined with KamLAND have measured the so- 
called solar parameters 1 Am^ = Am^ = 7t 3 -10~ 5 eV 2 and sin 2 (2# so i) = sin 2 (26>i 2 ) = 
0.8±g;l |Mal03l ICle98l ISAC02I IcTAT^I ICNOOOl ISKn2hl ISN()n2l IHuh02| . Future 
solar neutrino data as well as the forthcoming KamLAND results will undoubt- 
edly improve the solar neutrino parameters determination. A new middle baseline 
(20-70 km) reactor neutrino experiment could further improve Am^ or/and 6 so \ if 
necessary |Pet02l ISch03| . 

Atmospheric neutrino experiments such as Super-Kamiokande together with the 
K2K first long baseline accelerator neutrino experiment have measured the so-called 
atmospheric parameters Am% tm = |Am 32 | = 2+ ; 7 • 10~ 3 eV 2 and sin 2 (2# atm ) = 
sin 2 (2(9 32 ) = 1.0t£° | SK981 IMATMl IRonOll ISK02aj . Experimental errors will 
slowly decrease with additional K2K and Super-Kamiokande data, but a major 
improvement of the results is expected from the currently starting MINOS long 
baseline neutrino experiment |MIN01allMIN01b| . 

The third sector of the neutrino oscillation matrix is driven by the mixing angle # 13 , 
currently best constrained by the CHOOZ reactor neutrino experiment CH098 
ICHQ99I OKlOOl ICMO03| . CHOOZ provides the upper bound sin 2 (20 13 ) < 0.20 
(90 % C.L.), assuming Am 2 tm = 2.0 10~ 3 eV 2 (this upper limit is strongly corre- 
lated with the assumed value of Am^ tm .) A weaker upper bound, sin 2 (20i 3 ) < 0.4, 
has been obtained by the Palo- Verde experiment |PV01| . 



Concerning the determination of the PMNS mixing parameters, the measure- 
ment of the angle 613 is the next experimental step to accomplish. Knowing the 
value of 9is, or lowering the CHOOZ bound is already fundamental, in itself, in or- 
der to better understand the structure of the PMNS matrix. Both atmospheric and 



lr The intervals vary slightly in the different analyzes, we give here the values quoted in Hub02 . 
Furthermore, we assume here the normal neutrino mass hierarchy case. 
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solar mixing angles have been found to be maximal or large, thus the smallness of 
#13 remains a mystery. Moreover, any sub-leading three-neutrino oscillation effects, 
such as the solar-atmospheric driven oscillation interferences |Pet02l ISch03| or the 
CP-violation in the lepton sector, could only be observable for non- vanishing # 13 
values. 

Which sensitivity is then relevant for the forthcoming projects dedicated to 
#13 ? On the one hand, neutrino mass models predict sin 2 (2^13) values ranging 
from to 0.18 RWP04 . Any neutrino experiment with a sensitivity of a few per- 
cents, like Double-CHOOZ, has thus an important discovery potential. On the other 
hand, the neutrino mass models connect, in most cases, the CPS phase to the lepto- 
genesis mechanism |Buc04| . The search for CP violation effects in the lepton sector 
is thus of great interest since the leptogenesis mechanism is one of the best current 
explanation of the matter-antimatter baryon asymmetry observed in our Universe. 
The target sensitivity to achieve is thus strongly driven by the potential of future 
CPS neutrino appearance experiments. In the distant future, CP violation could 
be observed at neutrino factories if sin 2 (26*13) > 0.001. However, on a shorter time 
scale, a value of sin 2 (26*i 3 ) of a few percent might allow superbeam based experi- 
ments, possibly combined with a large reactor neutrino detector, to probe part of 
the 5 — 6>i3 parameter space |Hub02| . 

Although they are not designed to measure #13 , a marginal improvement of the 
CHOOZ constraint can be obtained with conventional neutrino beams. For instance, 
the MINOS experiment jMINOla] may achieve a sensitivity sin 2 (26*13) < 0.1, while 
the CNGS experiments, OPERA and ICARUS |OiMO02al ICNG02bl lHuE04| . may 
improve the CHOOZ bound down to sin 2 (26»i 3 ) < 0.14 and sin 2 (26>i 3 ) < 0.09, 
respectively, if no excess of u e is observed after five years of data taking 2 (Am 2 tm = 
2.01Q- 3 e V 2 , 90 % C.L. ). The quoted values reduce to 0.05 jMTNOIhj . 0.08 [QPE03 
ICNG02b| . 0.04 |ICA02j . respectively, by neglecting matter effects, CP-<5 p hase (se t 
to zero), and mass hierarchy induced correlations and degeneracies Min02 H ub02| . 

Concerning the future of neutrino physics, the next generation of accelerator 
neutrino experiments coupled with powerful neutrino beams (the so-called Super- 
beam long baseline experiment) are primarily dedicated to the determination of 
the PMNS mixing matrix elements #13 and CP-<5, as well as the precise mea- 
surement of the atmospheric mass splitting and mixing angle, and the identifi- 
cation of the neutrino mass hierarchy (the sign of Am 2 2 ). After five years of 
data taking, the T2K experiment aims to reach the sensitivity sin 2 (26*13) < 0-02 
(90 % C.L.) T2K02 T2K031 IHuE02] : a similar sensitivity is foreseen by the NuMI 
Off- Axis project 3 |NUM02| . 

The observation of a v e excess in an almost pure neutrino beam at any 
accelerator experiment would be major evidence for a non- vanishing #13. But un- 
fortunately, in addition to the statistical and systematic uncertainties, correlations 
and degeneracies between #13, 6> a tm, sgn(Am§ 1 ), and the CP-<5 phase degrade the 
knowledge of 6*i 3 |Min02l lHub02] . Even though appearance experiments seem to 
be the easiest way to measure very small mixing angles, as might be the case for 
#13, it is of great interest nevertheless to get additional information with another 
experimental method. 

A reactor neutrino experiment, like Double-CHOOZ, is able to measure #13 with 
an independent detection principle (inverse neutron beta decay), and thus different 

2 ICARUS and OPERA results could be combined, leading to a value very close to the ICARUS 
sensitivity (10 % improvement). 

3 This value takes into account, in a very conservative manner, all correlations and degeneracies. 
At a fixed 8 phase taken to be 0, the value quoted would be three times lower. 
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systematic uncertainties. Unlike appearance experiments, it does not suffer from 
parameter degeneracies induced by the CPS phase. In addition, thanks to the low 
V e energy (a few MeV) as well as the very short baselines (a few kilometers) the 
reactor measurement is not affected by matter effects. As a consequence reactors 
provide a clean information on sin 2 (26*13). Doublc-CHOOZ will use two identical 
detectors at ~150 m and 1.05 km from the CHOOZ-B nuclear power plant cores. 
The near detector is used to monitor both the reactor 77 e flux and energy spectrum, 
while the second detector is dedicated to the search for a deviation from the expected 
(1/distance) 2 behavior, tagging an oscillation effect. For Am 2 tm = 2.0 10~ 3 eV 2 wc 
expect a sensitivity of sin 2 (26*13) < 0.03 (90 % C.L.) after three years of data taking. 

In conclusion, due to the fundamental interest of #13 as well as the importance 
of its amplitude for the design of future neutrino experiments dedicated to CP- 
<5, independent ^-dedicated experiments are mandatory. To accomplish this goal, 
both reactor and accelerator programs should provide the required independent and 
complementary results (RWP04j . 



Chapter 2 

Searching for sin 2 (2^13) with 
reactors 

2.1 Neutrino oscillations 

Neutrino flavor transitions have been observed in atmospheric, solar, reactor and 
accelerator neutrino experiments. To explain these transitions, extensions to the 
minimal Standard Model of particle physics are required. The simplest and most 
widely accepted extension is to allow neutrinos to have masses and mixing, similar 
to the quark sector. The flavor transitions can then be explained by neutrino 
oscillations. 

2.1.1 Quark mixing 

The Wolfenstein parameterization of the CKM matrix PDGOO is based on the very 
small mixing between the quarks. The mixing matrix is almost the identity matrix 
with only small corrections for the off-diagonal elements. It uses the observed 
quark mixing angles hierarchy 1 to introduce an expansion parameter A describing 
the mixing between u and s quarks. This leads to the parameterization 

/ 1 - i A 2 A AX 3 (p-i v ) \ 

Vckm^ -A 1-±A 2 AX 2 +G(A 4 ), (2.1) 

\ AA 3 (1 -p + ir,) -AX 2 1 / 

where A corresponds to the Cabibbo angle sin^c ~ 0.22, and the other parameters 
are roughly A ~ 0.83, p ~ 0.23 and 77 ~ 0.36 PDGOO . The latter parameter 
describes CP violation in the quark sector; all such effects are proportional to |Jar85| 

J CP ~ -A 2 A 6 T] ~ -3 • 10~ 5 . (2.2) 
Therefore, CP violation in the quark sector is a small effect. 

2.1.2 Neutrino mixing 

The neutrino oscillation data can be described within a three neutrino mixing 
scheme, in which the flavor states v a (a = e, fx, r) are related to the mass states z/, 
(i = 1,2,3) through the PMNS (Pontecorvo-Maki-Nakagawa-Sakata) unitary lep- 
ton mixing matrix. 

1 e 12 - 0.1 > 023 ~ 0.01 > 013 ~ 0.001 
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It can be parameterized as Upmns = 






Cl3 si3e 
1 

-i 3 e tS c 13 J V 1 J V eI/3 

(2.3) 

C13C12 C13S12 

-C23S12 - Sl3S23Cl2e il5 C23C12 - Sl3S23Sl2e 115 

S23S12 - si 3 c 2 3Ci 2 e lS -S23C12 - si 3 C23Si2e l5 C13C23 / \ e l13 

where Cjj = cos#.y and = sin&y, 5 is a Dirac CP violating phase, a and /? 
are Majorana CP violating phases, not considered in the following. Up to now, 
the angles 9n and $23 are probed via the oscillations of solar/reactor and atmo- 
spheric neutrinos, while the angle #13 is mainly constrained by the CHOOZ reactor 
experiment; the Dirac phase 5 has not been constrained yet. 

The factorized form of this PMNS mixing matrix is often used to identify the 
mixing angles reported by the experiments 2 

#23 — #atni> #12 — #sol, and #13 = #chooz- (2-4) 
The relevant formula for the oscillation probabilities is 



P(„ a - vp ) = 8 a p - 2Re £ U m U* aj U Pj ^1 - exp , (2.5) 



Since the identification of the MSW-LMA mechanism as the solution of the solar 
neutrino anomaly |SK02bl ISN( )02l IK AM02j . we now know that the mass eigenstate 
with the larger electron neutrino component has the smaller mass (state 1). Solar 
neutrino oscillations occur then mainly together with the little heavier state 2: 

Am^! = ml - m\ = Am 2 ol > 0. (2.6) 

The large mass squared difference measured in the atmospheric sector is therefore 
the splitting between the mass eigenstate 3 and the more closely spaced 1 or 2. In 
addition, the CHOOZ reactor neutrino experiment shows that the mass eigenstate 
3 has only a very small electron neutrino component. In this description, the sign 
of the splitting between state 3 and states 1 and 2 is unknown; this leads to two 
possibilities of mass ordering: 

|Am 2 2 | = |m 2 -m 2 | = Am 2 tm . (2.7) 

Thus, one defines the normal hierarchy (NH) scenario m 3 > to 2 > mi, and the 
inverted hierarchy scenario (IH) m 2 > mj > m 3 . The determination of the sign of 
Am§ 2 is one of the next goals in neutrino oscillation physics. 



2.2 Measurement of sin 2 (2^i 3 ) with reactor v e 
2.2.1 Reactor V e flux 

The fissionable material in the CHOOZ pressurized water reactors (PWR) mainly 
consists of 235 U and 239 Pu, which undergo thermal neutron fission. The fresh fuel is 

2 Thanks to the smallness of f^ m 2 ao1 and sin 2 Schooz- 
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enriched to about 3.5 % in 235 U. Fast fission neutrons are moderated by light water 
pressurized to 150 bar. The dominant natural uranium isotope, 238 U, is fissile only 
for fast neutrons (threshold of 0.8 MeV) but it also generates fissile 239 Pu by thermal 
neutron capture, 

n + 238 U ^ 239 U ^ 239 Np ^ 239 Pu (T 1/2 = 24, 100 y). (2.8) 

The 2 Pu isotope is produced in a manner similar to 239 Pu 

n + 239 Pu ^ 240 Pu + n ^ 241 Pu (T 1/2 = 14.4 y). (2.9) 

As the reactor operates, the concentration of 235 U decreases, while that of 239 Pu 
and 241 Pu increases. After about one year, the reactor is stopped and one third 
of the fuel elements are replaced. Typical numbers for an annual cycle are given 
in Table |2~T1 Due to the threshold of the detection reaction at 1.8 MeV only the 





Mean energy per fission 


Refueling 


cycle 




(MeV) 


beginning 


end 


235TJ 


201.7 ± 0.6 


60.5 % 


45.0 % 


238tj 


205.0 ± 0.9 


7.7 % 


8.3 % 


239p u 


210.0 ± 0.9 


27.2 % 


38.8 % 


241p u 


212.4 ± 1.0 


4.6 % 


7.9 % 



Table 2.1: Typical fuel composition for an annual cycle of a PWR power station, 
for the four main isotopes, normalized to 100 %. There are also other isotopes, not 
included here, which contribute for a few percents. 

most energetic antineutrinos are detected; they correspond to the decay of fission 
products with the highest Q-values and hence to the shortest lived. The detected 
antineutrinos thus closely follow changes in power. In particular spent fuel elements 
which are kept on site out of the core contain only long lived emitters with a low Q- 
value; their contribution to the detected v c signal is negligible. Measurements of the 
neutrino rate per fission have been performed for 235 U, 239 Pu and 241 Pu by Borovoi 
et al. [Bor83| and Schreckenbach et al. Sch85, Hah89j. The latter measurement 
includes the shape of the energy spectrum, with a 2 % bin-to-bin accuracy and an 
overall normalization error of 2.8 %. The measurement performed by |Sch85| can 
be compared with several computations and is found to be in good agreement with 
that of |Kla82l|Vog81| . We will therefore use this computation for the 238 U neutrino 
rate, which has never been measured. The 238 U contribution to the total number 
of fissions is ^10 %, and is therefore not a major source of error. The antincutrino 
spectra of the four dominant fissioning isotopes are shown in Figure lTTI During the 
cycle, the contributions of the different fissile isotopes to energy production evolve. 
For fresh fuel, 235 U fissions dominate, whereas 238 U fissions amount for a few times 
less. Quickly after the beginning of the cycle 239 Pu gives an important contribution 
(see Figure |2~2|) . 

2.2.2 77 e detection principle 

Reactor antineutrinos are detected through their interaction by inverse neutron 
decay (threshold of 1.806 MeV) 

v e +P -* e+ + n . (2-10) 
The cross section for inverse /3-decay has approximately the form 

2tt 2 ^ 3 

a(E e+ ) ~ ^—p e+ E e+ , (2.11) 
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Figure 2.1: v c spectra of the four dominant isotopes with their experimental error 
bars ( 238 U spectrum has not been measured but calculated). 



where p c + and E c + are the momentum and the energy of the positron 3 , r n is the 
lifetime of a free neutron and / is the free neutron decay phase space factor. As an 
approximation, we use an averaged fuel composition typical during a reactor cycle 
corresponding to 235 U (55.6 %), 239 Pu (32.6 %), 238 U (7.1 %) and 241 Pu (4.7 %). 
The mean energy release per fission W is then 203.87 MeV and the energy weighted 
cross section amounts to 

< c >fission= 5.825 • 10~ 43 cm 2 per fission . ( 2 -12) 

The reactor power Pth is related to the number of fissions per second Nf by 

N f = 6.241 ■ lO^sec" 1 • {P th [MW])/(W[MeV]) . (2.13) 

The event rate at a distance L from the source, assuming no oscillations, is thus 

i? L - N t ■ < a > fission • n p ■ \/{AttL 2 ) , (2.14) 

where n p is the number of protons in the target. For the purpose of simple scaling, a 
reactor with a power of 1 GWtd induces a rate of ^450 events per year in a detector 
containing 10 29 protons, at a distance of 1 km. 

Experimentally one takes advantage of the coincidence signal of the prompt 
positron followed in space and time by the delayed neutron capture. This very clear 
signature allows to strongly reject the accidental backgrounds. The energy of the 
incident antineutrino is then related to the energy of the positron by the relation 

£k - E c+ + (m„ - m p ) + 0(EhJm n ) . (2.15) 

Experimentally, the visible energy seen in the detector is given by E v i s = E c + + 
511 keV, where the additional 511 keV come from the annihilation of the positron 
with an electron when it stops in the matter. 

3 E e + is the sum of the rest mass and kinetic energy of the positron. 
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Days 



Figure 2.2: Percentage of fissions of the four dominant fissile isotopes during 300 
days of a typical fuel cycle. 



2.2.3 i> e oscillations 

Reactor neutrino experiments measure the survival probability Py a —iu B of the elec- 
tron antineutrinos emitted from the nuclear power plant 4 . This survival probability 
does not depend on the 5- CP phase. Furthermore, because of the low energy as well 
as the short baseline considered, matter effects are negligible |Min02| . Assuming a 
"normal" mass hierarchy scenario, m\ < m 2 < m^, the T> c survival probability can 
be written |Bil01| IPetOlj 

iW. = l-2sin 2 13 cos 2 13 sin 2 (^^) (2.16) 

l/i q , o / Am 2 ,! L 

- -cos 4 0i 3 sin 2 (26>i 2 )sin 2 ' 



2 v 7 V 4E 
+ 2 sin 2 fl 13 cos 2 0i3 sin 2 12 ( cos ( — ^ -cos ^| 

The first two terms in Ea. l2.2.3*l contain respectively the atmospheric driven (Am 31 = 
Am 2 tm ) and solar driven (Am 21 = Am^, 6yi ~ $soi) contributions, while the third 
term, absent from an?/ two-neutrino mixing model, is an interference between solar 
and atmospheric driven oscillations whose amplitude is a function of 613. Thus, up 
to second order in sin 2#i 3 and a — f"" 3 °' the survival probability can be expressed 
as 

Pv B ^p e ^ 1 - sin 2 26»i3 sin 2 (Am^L/4£;) + a 2 (Am 2 31 L/AE) 2 cos 4 6 13 sin 2 26 12 , 

(2.17) 



'The low neutrino energy (a few MeV) does not allow any appearance measurement. 
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where the third term on the right side can safely be neglected given the current 
range (90 % error intervals) of mixing parameters found in neutrino oscillation 
experiments 5 |SK02al EKMa\ : 

( A TOatm)sK-I = 

(sin 2 20 23 ) S k-i = 

( ATO atm)sK-L/E = 

(sin 2 26> 23 )sk-l/e = 
sin 2 (20 12 ) = 

Reactor experiments thus provide a clean measurement of the mixing angle #13, 
free from any contamination coming from matter effects and other parameter corre- 
lations or degeneracies [ MinQ2l fHub02] . Therefore they are exclusively dominated 
by statistical and systematic errors. 



2.0± 7 • 10" 3 eV 2 

1+0 
x -o.i 

2.4±° j • 10~ 3 eV 2 

1+0 
i -o.i 

7.0t 2 3 ■ 10" 5 eV 2 



2.3 Complementarity with Superbeam experiments 

A very detailed comparison of reactor antincutrino experiments with superbeams 
is described in |Min02l [Hub02 . Forthcoming accelerator neutrino experiments, or 
superbeams, will search for a v e appearance signal. The appearance probability 
Py -^Ve with terms up to second order, e.g., proportional to sin 2 2#i3, sin 20i3 • a, 
and a 2 , can be written as: 

P v ^ u , sin 2 20i3 sin 2 23 sin 2 {Am 2 31 L/AE) 

=F a sin 2#i 3 sin 8 sin 20 i2 sin 20 23 

(Am^L/4£) sin 2 (Am^L/4B) 
— a sin 20i3 cos 5 sin 2(9i 2 sin 20 2 3 

{kmlxL/AE) cos (Am^I/4B) sin (Am^I/4B) 
+ a 2 cos 2 23 sin 2 20i 2 (Am| 1 i/4£;) 2 , (2.18) 

where the sign of the second term refers to neutrinos (minus) or antineutrinos (plus). 
From Equation 12.181 one sees that superbeams suffer from parameter correlations 
and degeneracies coming from the different combinations of parameters. Many of the 
degeneracy problems originate in the summation of the four terms in Eauation l2.18l 
since changes in one parameter value often can be compensated by adjusting another 
one in a different term. This leads to the (5, $13) |Bur01) . sgn(Am| 1 ) .MinOTj , and 
(023, 7r/2 — 23 ) |Fog98| degeneracies, e.g. an overall "eight-fold" degeneracy [Bar02 
IMinOlj . Table O 

summarizes the sensitivity of accelerator and Doublc-CHOOZ 

experiments. 



5 Two different best fit values for the atmospheric mass splitting have been released by the 
Super-Kamiokande collaboration, based on two different analyzes of the same data. 
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Chooz Beams Double-CHOOZ 


T2K 


sm 2 (26»i 3 ) sensitivity limit (90 % CL) 


sin 2 (26ii3) 0.2 0.061 0.032 


0.023 


sin 2 (20i 3 )eff 0.2 0.026 0.032 


0.006 


Measurements for large sm J (29 13 ) = 0.1 (90 % CL) 


c;„ 2 roa ^ n 1+0.104 n 1+0.034 


n 1+0.067 
u ' 1 -0.034 



Table 2.2: Comparison of the sensitivity of reactor and accelerator based future 
neutrino experiments. The results of the table have been extracted from |Hub04j . 
"Beams" is the combination of the forthcoming MINOS, ICARUS, and OPERA 
experiments. Results for accelerator experiments are given for five years of data 
taking. Results for Double-CHOOZ are given for three years of operation. The line 
starting by "sin 2 (26*13)" provides the results of the computation taking into account 
all correlation and degeneracy effects, while the line starting by sin 2 (26 , i3) e ff give 
the results of a similar computation performed after "switching off" those effects. 
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Chapter 3 



Overview of the 
Double-CHOOZ experiment 

This section is an overview of the Double-CHOOZ experiment. The Double-CHOOZ 
technology of reactor neutrino detection is based on experience obtained in numer- 
ous experiments: Goesgen GOE86:, Bugey BUG96 (at short distances), CHOOZ 
|GH098l \TWm IOHOOOI ICTinOSl . Palo Verde |FVoTl (at km scale distance) 
CTF |CTF98| . Borexino |Sch99| . Kamland |KAM02| (distances of a few hundred 
km). Therefore, no long term R&D program has to be conducted prior to designing 
and building the new detector. Nevertheless, in order to be a precision experiment, 
the Double-CHOOZ design has to be improved with respect to CHOOZ. The liquid 
scintillator is described in Chapter the calibration in Chapter El and the back- 
grounds in Chapter The systematic error handling is presented in Chapter [SJ To 
conclude, the sensitivity of Double-CHOOZ, taking into account the overall set of 
systematic errors, is presented in Chapter |2l 

3.1 The F e source 

To fulfill the aim of the Double-CHOOZ experiment, precise knowledge of the an- 
tineutrino emission in the nuclear core is not crucial thanks to the choice of compar- 
ing two similar detectors at different distances, where the near detector measures 
the flux without v c losses consequent to oscillations. Nevertheless this information 
is available and will be used for cross checks and other studies of interest (see Ap- 
pendix^). More details on the V e source will be necessary to do those specific 
studies with the near detector, such as the v c spectrum and flux expected for a 
given fuel composition and burn up. 

3.1.1 The CHOOZ nuclear reactors 

The antineutrinos used in the experiment are those produced by the pair of reactors 
located at the CHOOZ-B nuclear power station operated by the French company 
Electricite de France (EDF) in partnership with the Belgian utilities Electrabel 
S.A./N.V. and Societe Publique d'Electricite. They are located in the Ardennes 
region, northeast of France, very close to the Belgian border, in a loop of the Meuse 
river (see Figures mi and ET^ . At the CHOOZ site, there are two nuclear reactors, 
both are of the most recent N4 type (4 steam generators) with a thermal power 
of 4.27 GWth, and recently upgraded from 1.45 GW C to 1.5 GW G . These reactors 
are of the Pressurized Water Reactor type (PWR) and are fed with UOx type fuel. 
They are the most powerful reactor type in operation in the world. One unusual 
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characteristic of the N4 reactors is their ability to vary their output from 30 % 
to 95 % of full power in less than 30 minutes, using the so-called gray control 
rods in the reactor core. These rods are referred to as gray because they absorb 
fewer free neutrons than conventional black rods. One advantage is a bigger thermal 
homogeneity. 205 fuel assemblies are contained within each reactor core. The entire 
reactor vessel is a cylinder of 13.65 meters high and 4.65 meters in diameter. The 
first reactor started operating at full power in May 1997, and the second one in 
September of the same year. 

3.2 Detector site 

The Double-CHOOZ experiment will run two almost identical detectors of medium 
size, containing 12.7 cubic meters of liquid scintillator target doped with 0.1 % of 
Gadolinium (see Chapter 0). The neutrino laboratory of the first CHOOZ experi- 
ment, located 1.05 km from the two cores of the CHOOZ nuclear plant will be used 
again (see Figure EOjl . This is the main advantage of this site compared with other 
French locations. We label this site the far detector site or CHOOZ-far. A sketch 




Figure 3.1: Overview of the experiment site. 



of the CHOOZ-far detector is shown in FigureEH The CHOOZ-far site is shielded 
by about 300 m.w.e. of 2.8 g/cm 3 rocks. Cosmic ray measurements were made with 
Resistive Plate Chambers and compared with the expected angular distributions. 
A geological study revealed the existence of several very high density rock layers 
(3.1 g/cm 3 whose positions and orientations were in agreement with the cosmic ray 
measurements [CHO03] ). It is intended to start taking data at CHOOZ-far at the 
beginning of the year 2007. 

In order to cancel the systematic errors originating from the nuclear reactors 
(lack of knowledge of the 77 c flux and spectrum), as well as to reduce the set of sys- 
tematic errors related to the detector and to the event selection procedure, a second 
detector will be installed close to the nuclear cores. We label this detector site the 
near site or CHOOZ-near. Since no natural hills or underground cavity already 
exists at this location, an artificial overburden of a few tens of meters height has to 
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Figure 3.2: Map of the experiment site. The two cores are separated by a distance 
of 100 meters. The far detector site is located at 1.0 and 1.1 km from the two cores. 



be built. The required overburden ranges from 53 to 80 m.w.e. depending on the 
near detector location, between 100 and 200 meters away from the cores (see Ta- 
ble l3.1(l . A sketch of this detector is shown in Figure l3~5l After first discussions, this 



Distance 


Minimal overburden 


Required overburden 




(m.w.e.) 


(m.w.e.) 


100 


45 


53 


150 


55 


65 


200 


67.5 


80 



Table 3.1: Overburden required for the near detector. The second column is the es- 
timation of the minimal overburden required for the experiment. The third column 
is minimal overburden added to a safety margin. 

construction has been considered as technically possible by the power plant com- 
pany authorities. An initial study has been commissioned by the French electricity 
power company EDF to determine the best combination of location-overburden and 
to optimize the cost of the project. Plan is to start to take data at CHOOZ-near 
at the beginning of the year 2008. 



3.3 Detector design 

The detector design foreseen is an evolution of the detector of the first experi- 
ment |(JH()03| . To improve the sensitivity of Double-CHOOZ with respect to 
CHOOZ it is planned to increase statistics and to reduce the systematic errors 
and backgrounds. 

In order to increase the exposure to 60,000 V e events at CHOOZ-far (statistical 
error of 0.4 %) it is planned to use a target cylinder of 120 cm radius and 280 cm 
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Figure 3.3: Picture of the CHOOZ-far detector site taken in September 2003. The 
original CHOOZ laboratory hall constructed by EDF, located close the the old 
CHOOZ-A underground power plant, is still in perfect condition and could be re- 
used without additional civil engineering construction. 



height, providing a volume of 12.7 m 3 , ~2.5 larger than in CHOOZ. In addition, 
the data taking period will be extended to at least three years, and the overall data 
taking efficiency will be improved. The global load factor of the reactor, i.e. the 
average reactor efficiency, is about 80 %, whereas it was significantly lower for the 
CHOOZ experiment performed during the power plant commissioning. In addition, 
the detector efficiency will be slightly improved. The background level at CHOOZ- 
far will be decreased to have a signal to noise ratio over 100 (about 25 in CHOOZ). 

The near and far detectors will be identical inside the PMTs supporting struc- 
ture. This will allow a relative normalization systematic error of ^0.6 % (see 
Chapter |HJ. However, due to the different overburdens (60-80 to 300 m.w.e.), 
the outer shielding will not be identical since the cosmic ray background varies 
between CHOOZ-near and CHOOZ-far. The overburden of the near detector has 
been chosen in order to keep the signal to background ratio above 100. Under this 
condition, even a knowledge of the backgrounds within a factor two keeps the associ- 
ated systematic error well below the percent (assuming that its energy distribution 
is known). 

The detector design has been intensively studied and tested with Monte-Carlo 
simulations, using two different softwares derived from the simulation of the CHOOZ 
and Borexino experiments (see Chapter |2J . In order to increase the width of the 
liquid buffers protecting the v c target, the 1 meter low radioactivity sand shielding 
of CHOOZ will be replaced by a 15 cm metal shielding, steel or iron (this is used 
to reduce the external gamma rays coming from the rock.) This will increase the 
size of the liquid active buffer and will thus improve the rejection of muon induced 
backgrounds (see Chapter d) . Starting from the center of the target the detector 
elements are as follows (see Figures 13.41 and l3~5*f 

• F e target 
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Figure 3.4: The new CHOOZ-far detector, at the former CHOOZ underground site. 
The detector is located in the tank used for the CHOOZ experiment (7 meters high 
and 7 meters in diameter) that is still available. About 12.7 m 3 of a dodecane+PXE 
based liquid scintillator doped with gadolinium is contained in a transparent acrylic 
cylinder surrounded by the 7-catcher region and the buffer. The design goal is 
to achieve a light yield of about 200 pe/MeV (see Chapter 0J) which requires an 
optical coverage of about 15 %, provided by the surrounding PMTs. The PMTs 
are mounted on a cylindrical structure which separates optically the outer part of 
the detector, which is used as a muon veto. 



A 120 cm radius, 280 cm height, 6-10 mm width acrylic cylinder, filled with 
0.1 % Gd loaded liquid scintillator target (see Chapter |5J). 

• 7-catcher 

A 60 cm buffer of non-loaded liquid scintillator with the same optical proper- 
ties as the v e target (light yield, attenuation length), in order to get the full 
positron energy as well as most of the neutron energy released after neutron 
capture. 

• Buffer 

A 95 cm buffer of non scintillating liquid, to decrease the level of accidental 
background (mainly the contribution from photomultiplier tubes radioactiv- 

ity). 

• PMT supporting structure 

• Veto system 

A 60 cm veto region filled with liquid scintillator for the far detector, and a 
slightly larger one (about 100 cm) for the near detector. 
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Figure 3.5: The CHOOZ-near detector at the new underground site, close to the 
reactor cores. This detector is identical to the CHOOZ-far detector up to and 
including the PMT surface. The veto region will be enlarged to better reject the 
cosmic muon induced backgrounds (see Chapter 0. 



Compared to previous scintillator 77 G detectors, the Double-CHOOZ experiment will 
use cylindrical targets; Monte-Carlo simulation shows that the spatial reconstruc- 
tion in a cylinder is suitable for the experiment. A spherical configuration gives 
slightly better results, however. Each parameter of the detector is being studied by 
Monte-Carlo simulation in order to define the tolerance on the differences between 
the two detectors (see Chapter 0J . The inner volume dimensions as well as the 
shape of the target vessels are still preliminary, within a few tens of percents, and 
could change prior to the publication of the proposal. 

3.3.1 Experimental errors and backgrounds 

In the first CHOOZ experiment, the total systematic error amounted to 2.7 % CHO03 . 
Table 13.21 summarizes the control of the systematic uncertainties that had been 
achieved in the first CHOOZ experiment as well as the goal of Double-CHOOZ. 
The main uncertainties at CHOOZ came from the 2 % only knowledge of the an- 
tineutrino flux coming from the reactor. This systematic error vanishes by adding 
a near detector to monitor the power plant antineutrino flux and energy spectrum. 
A complete description of the systematic uncertainties is given in Chapter |H1 The 
main challenge of the Double-CHOOZ experiment is to decrease the overall system- 
atic error from 2.7 % to 0.6 %. The strategy is to improve the detector design, to 
rely on the comparison of the two detectors, and to reduce the number of analysis 
cuts. The non-scintillating buffer will reduce the singles rates in each detector by 
two orders of magnitude with respect to CHOOZ. This allow to lower the positron 
threshold down to ^500 keV, well below the 1.022 MeV physical threshold of the 
inverse beta decay reaction. A very low threshold has three advantages: 

• The systematic error due to this threshold is suppressed. It was one of the 
largest source of systematic error, 0.8 % in CHOOZ CHO03 . 

• The background below the physical 1 MeV threshold can be measured. 
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• The onset of the positron spectrum provides an additional calibration point 
between the near and far detectors. 

This reduction of the singles events relaxes or even suppresses the localization cuts, 
such as the distance of an event to the PMT surface and the distance between 
the positron and the neutron. These cuts, used in CHOOZ jCHO03| . are difficult 
to calibrate and have to be avoided or relaxed in Double-CHOOZ. The remaining 
event selection cuts will have to be calibrated between the two detectors with a 
very high precision. Most important will be the calibration of the energy selection 
of the delayed neutron after its capture on a Gd nucleus (with a mean energy 
release of 8 MeV gammas). The requirement is ~100 keV on the precision of 
this cut between both detectors, which is feasible with standard techniques using 
radioactive sources (energy calibration) and lasers (optical calibrations) at different 
positions throughout the detector active volume (see Chapter 0) . The sensitivity 
of a reactor experiment of Double-CHOOZ scale (~300 GWth .ton. year) is mostly 
given by the total number of events detected in the far detector. The requirement 
on the positron energy scale is then less stringent since the weight of the spectrum 
distortion is low in the analysis. (This is being studied by simulation.) A detailed 



CHOOZ Double-CHOOZ 



Reactor cross section 


1.9 % 




Number of protons 


0.8 % 


0.2 % 


Detector efficiency 


1.5 % 


0.5 % 


Reactor power 


0.7% 




Energy per fission 


0.6 % 





Table 3.2: Summary of the systematic errors in CHOOZ and Double-CHOOZ (goal). 
The first line, "Reactor cross section" , accounts for the uncertainties on the neutrino 
flux as well as the inverse neutron beta decay cross section. A two F e detectors 
concept makes the experiment largely insensitive to the "Reactor cross section" and 
the reactor power uncertainties. The number of protons in the first acrylic vessel 
targets as well as the detection efficiencies have then to be calibrated between the 
two detectors, but only in a relative sense. 



background study is presented in Chapter [7| In CHOOZ the dominant correlated 
proton recoil background was measured to be about one event per day |CHO03j . At 
CHOOZ-far the active buffer will be increased, with a solid angle for the background 
being almost unchanged. This together with a signal increased by about a factor 
of 3 will fulfill the requirement of a signal to noise ratio greater than 100. At 
CHOOZ-near, due to the shallow depth between 60 and 80 m.w.e, the cosmic ray 
background will be more important. If, for instance, the CHOOZ-near detector is 
located 150 meters from the nuclear cores, the signal will be a few thousand events 
per day, while the muon rate is expected to be a factor of ten less. A dead time 
of about 500 /xsec will be applied to each muon 1 , leading to a global dead time of 
about 30 %. A few tens of recoil proton events per day, mimicking the V c signal, 
are expected while the estimate of the muon induced cosmogenic events ( 9 Li and 
8 He) is less than twenty per day with a large uncertainty (this last point is being 
carefully studied). This fulfills the requirement of a signal to noise ratio greater 
than 100 at CHOOZ-near. 



This is a conservative number that could be reduced. 
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Near and Far Detectors Spectra 
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Figure 3.6: Positron spectrum (visible energy, MeV) simulated for the CHOOZ-near 
and far detectors 



3.3.2 Sensitivity 

A detailed study of the Double-CHOOZ sensitivity is presented in Chapter^ From 
the simulations, we expect a sensitivity of sin 2 (26*13) < 0.03 at 90 % C.L. for 
Am 2 tm = 2.0 10~ 3 eV 2 (best fit value of Super-Kamiokande |SK02ap . after three 
years of operation. According to the latest Super-Kamiokande L/E analysis the best 
mass splitting is Am 2 tm = 2.4 10~ 3 eV 2 |SK04a| . The Double-CHOOZ sensitivity 
would then be sin 2 (26*13) < 0.025 2 . A study of the evolution of the sensitivity with 
respect to the luminosity is presented in Figure 1X71 Hub04 . A sensitivity of ~0.05 
is reachable within the first year of operation with 2 detectors. These estimates are 
based on the assumptions that the relative normalization error between the near 
and far detectors could be kept at 0.6 %, and that the backgrounds at both sites 
amount to 1 % of the V c signal (we assume those backgrounds to be known within 
a factor of two). 

The effect of v c oscillations on the positron spectrum is displayed in Figure l3~%l 
for different values of Am 2 tm and sin 2 (26*13). For Am 2 tm > 2.0 10~ 3 eV 2 a sig- 
nificant shape distortion is expected at the onset of the energy spectrum. As- 
suming sin 2 (26 ) i3) = 0.15, the ratio of the near and far detector spectrum is pre- 
sented in Figure 13.91 with the expected statistical error bars (1 a) after three 
years of data taking. It is worth mentioning that the 1.05 km average base- 
line at CHOOZ is not optimal (the optimal distance should be roughly 1.5 km) 
compared to the first oscillation maximum if the atmospheric mass splitting is 
Am 2 tm = 2.0 10 -3 eV 2 . Nevertheless, a new Super-Kamiokande analysis of the 
data indicates 1.9 10" 3 eV 2 < Am 2 tm < 3.0 10" 3 eV 2 (90 % C.L.), with a best fit 
at Am 2 tm = 2.4 10 -3 eV 2 }SK04a} ■ A shorter baseline is compensated by higher 



2 The sensitivity of the other experiments is as well better for higher Am 2 
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Far detector at L=1.05 km 
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Figure 3.7: Luminosity scaling of the Double-CHOOZ sin 2 (2^13) sensitivity at the 
90 % C.L.. Here, Aml tm = 2.0 10~ 3 eV 2 is assumed to be known within 5 %. The 
relative normalisation error between the two detectors is taken to be (0.2 %)0.6 % 
for the light (dark) shaded regions. Correlated backgrounds with known shapes 
account for 1 % and are supposed to be known within 50 %. A 0.5 % "Flat" bin-to- 
bin uncorrelated background component as been accounted as well (known within 
50 %). A luminosity of 300 GWth • ton • year (left vertical line) correspond approxi- 
mately to the setup of the Double-CHOOZ experiment as described in this Letter of 
Intent (sin 2 (20 13 ) < 0.03 within 3 years of data taking). However, a luminosity of 
8000 GWth • ton • year (right vertical line) would describe a ~300 tons next detector 
generation at Chooz |Hub04j 



statistics for a fixed size detector. However, a value of Am 2 tm < 1.5 10 3 eV 2 would 
restrict the absolute potential of the Double-CHOOZ experiment (see Chapter |5J|. 
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Figure 3.8: Ratio of expected number of V e events in the far detector with respect to 
the no oscillations scenario, after 3 years data taking, for different values of Am 2 tm 
and sin 2 (2^13). 
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Figure 3.9: Ratio of observed number of events in the far detector with respect 
to the no oscillations scenario, after 3 years data taking, for sin 2 (26*13) = 0.15 and 



different values of Ar 



The error bars plotted here are only statistical (ler). The 



positron spectrum shape is also displayed in the background. The potential of the 
experiment to exclude sin 2 (2f5 13 ) = may be seen as a deviation from unity in the 
ratio. Note that in some cases spectral information may be important. The largest 
spectrum deviation effect is located at the onset of the spectrum, below 4 MeV. 



Chapter 4 

Detector design and 
simulation 

In this section we describe the detector design envisaged for the Double-CHOOZ 
experiment. Although the generic design is almost complete, some specific technical 
solutions are still preliminary and could evolve prior to publication of the proposal. 

4.1 Detector design 

Detector dimensions are shown in Figure |4~T1 

4.1.1 The V e target acrylic vessel (12.67 m 3 ) 

The neutrino target is a 120 cm radius 280 cm height transparent acrylic cylinder, 
filled with 0.1 % Gd loaded liquid scintillator (see Chapter[SJl. Wall thicknesses (un- 
der study) range from 6 to 10 mm. The inner acrylic vessel is depicted in Figure FOl 
(left). Since the relative volume between the two inner acrylic vessels has to be 
controlled at a very accurate level (0.2 %), we plan to build both acrylic vessels at 
the manufacturer site and to move them as single units into the detector sites. This 
is possible for the far site, thanks to the size of the underground tunnel. The near 
detector site will be designed in order to allow this operation. With this strategy, no 
acrylic welding or gluing has to be done on site, thus reducing the uncontrolled dif- 
ferences between the two envelopes. Furthermore, a very precise calibration of both 
inner vessels is foreseen at the manufacturer (filling test). Current R&D focuses on 
the chemical compatibility between acrylic and liquid scintillator (see Chapter [SJ. 
Preliminary stress calculations have been done for this purpose (see Figure EH?)) . 



4.1.2 7-catcher acrylic vessel (28.1 m 3 ) 

The 7-catcher is a 180 cm radius and 400 cm height acrylic cylinder filled with 
non-loaded liquid scintillator, which has the same optical properties as the V e tar- 
get (light yield, attenuation length). Unlike the inner envelope, this second acrylic 
vessel will have to be partially assembled on site. Nevertheless, the shape and di- 
mensions between far and near 7-catcher are less critical than for the inner vessels. 
Therefore, small differences between the near and far 7-catcher acrylic vessels could 
be tolerated (a Monte-Carlo study is being done to provide the construction toler- 
ance). 
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Figure 4.1: Dimensions of the CHOOZ-far detector (in cm). Starting from the 
center we have: the neutrino target region composed of Gd doped liquid scin- 
tillator (12.7 m 3 ), the 7-catcher region composed of unloaded liquid scintillator 
(28.1 m 3 ), the non scintillating buffer region (100.0 m 3 ), and the veto (110.0 m 3 ). 
The CHOOZ-near detector is identical up to and including the PMT support struc- 
ture; however, its external muon veto is slightly larger to better reject the cosmic 
muon induced backgrounds. The exact PMT positioning has not been chosen yet. 



This scintillating buffer around the target is necessary to: 

Measure the gammas from the neutron capture on Gd. 

The total released energy is 8 MeV, with a mean gamma multiplicity of 3 to 4. 
But there are also some 8 MeV single gammas. The buffer must be thick enough 
to reduce the gamma escape out of the sensitive volume, i.e. the target and the 
7-catcher. This escape creates a tail below the 8 MeV peak. Since we must apply 
an energy cut to define the neutron capture on Gd, the tail of the energy spectrum 
has to be small enough to keep the systematic error negligible (if there is an energy 
scale mismatch between both detectors). Monte-Carlo simulations with a 60 cm 
buffer and a 100(150) keV energy error gives 0.2(0.3) % difference in the neutron 
counting, which is tolerable. 

Measure the positron annihilation. 

To have a clean threshold at 500 keV, it is mandatory to have very few events with 
an energy below 1 MeV. From the simulation, a thickness of 35 cm is adequate. 

Reject the background. 

This is the most demanding constraint. One of the most severe background in 
Double-CHOOZ comes from very fast neutrons, created by muons crossing rocks 



4.1. DETECTOR DESIGN 



35 




R= 180 

Figure 4.2: The two acrylic vessels containing the Gd doped and unloaded scintil- 
lators. The lines drawn on the cylinders show the preliminary positioning of the 
welding joint between the acrylic pieces. The inner envelopes will be constructed 
at the manufacturer and transported as single units to the detector sites while the 
outer envelopes will have to be assembled on sites. 




Figure 4.3: Preliminary evaluation of the stress applied on an empty acrylic cylin- 
der suspended with three kevlar ropes (set at 120 degrees from each other). The 
maximum stress has been estimated at 12 MPa, while acrylic supports a maximum 
of ^24 MPa in the elastic regime. 
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near the detector (see Chapter [7J| . To be able to reach the target traveling through 
the 2 m buffers, these neutrons must have an energy greater than 20 MeV. So 
when arriving at the scintillating buffer, they often deposit more than 8 MeV in 
the sensitive volume. This provides a useful rejection, by a factor of ~2. In the 
simulation, this rejection was seen to be stable for large buffer thickness, and to 
decrease when this thickness is reduced below 60 cm. Another advantage of this 
thickness is to allow to scale the result of the first experiment, since the sensitive 
volume around the target will be the same in both experiments (the veto volume 
was not sensitive to low energy events in the first experiment). 

4.1.3 Non scintillating buffer (100 m 3 ) 

The non-scintillating region aims to decrease the level of the accidental backgrounds, 
mainly due to the contribution from the photomultiplier tubes (see Chapter [JJ. To 
define the size of this region, we have to consider the following constraints: 

1. The fast neutron background implies to keep the distance from the rock to the 
neutrino target at least as it was for the CHOOZ experiment case 1 . Scaling 
from the CHOOZ experiment, we thus need at least 215 cm of liquid from the 
rock to the target. 

2. The size of the target has been chosen to be 120 cm to decrease the statistical 
error down to 0.4 %, after three years of operation. 

From those constraints, the total thickness of the veto and the non-scintillating 
buffer has to be smaller than 155 cm. Accounting for the size of the laboratory 
(mechanical constraints) and the requirement to have an efficient veto, we choose 
the thickness of the veto to be around 60 cm. From those considerations, the 
non-scintillating buffer region reduces to 95 cm. The simulation shows that this 
configuration fulfills our requirement on the accidental background level tolerated 
(which is mainly driven by PMT radioactivity). 



4.1.4 PMTs and PMT support structure 

The PMT support structure is a 275 cm radius and 590 cm height cylinder (material 
under study) filled with the same liquid as the 7-catcher, mixed with a quencher 
(DMP for instance). 

From the simulation, 500 PMTs of 8" are necessary to achieve a photoelectron 
yield of ^200 photoelectrons per MeV. Another possibility would be to use a smaller 
number of larger PMTs, 10" , 12" or 13" for instance. The reference PMT is the pho- 
tomultiplier 9354KB of ETL |ETT| . The glass used has a very low activity (60 ppm 
in K, 30 ppb in Th and U), and the quantum efficiency peaks at about 28 % at 
430 nm. For those PMTs, the peak-to- valley ratio of the single photoelectron signal 
is typically 2 (1.5 guaranteed by the manufacturer). Since we expect 600 photo- 
electrons for a medium energy signal of 3 MeV (visible energy), there will always 
be an important fraction of the PMTs working in the single photoelectron regime. 
The electronics gain is in the 10 6 — 10~ 7 range, hence some additional amplification 
is required in the front-end electronics system to obtain a good single photoelec- 
tron peak definition (additional dynodes could also be a way to increase the gain). 
Photonis as well as Hamamatsu PMTs are under study. The final photomultiplier 
choice will be made in 2004, during the design phase. 

1 The target vessel is seen from outside of the detector under a similar solid angle in both 
experiments. 
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Figure 4.4: Surface of PMTs mounted on the support structure of the detector as 
described in the GEANT4 simulation. About 500 PMTs are displayed. 



4.1.5 Veto (110 m 3 ) 

The external veto is contained in a steel cylinder of 350 cm radius and 710 cm 
height. The veto thickness is 60 cm for the far detector. It can be enlarged for 
the near detector, to better reject the cosmic muon induced backgrounds, since the 
laboratory has to be build. This tank is shielded by 15 cm of steel in order to reduce 
the external backgrounds. 

4.2 Fiducial volume 

4.2.1 Definition of the fiducial volume 

A neutrino interaction in this detector will be tagged by the neutron capture on 
gadolinium (as was the case in the first CHOOZ experiment CHO03 ). This is 
the main advantage of using a gadolinium loaded scintillator. However, there is an 
additional effect to consider, the spill in/out, that leads to a compensation between 
two kinds of v e interactions: 

• The v c interacts in the inner acrylic target, near the vessel, but the neutron 
escapes the target, and is captured on hydrogen in the 7-catcher. In that 
case, there is no Gd capture to characterize the neutrino interaction, and this 
is thus not selected as a neutrino event. 

• The T> c interacts in the 7-catcher, not too far from the target, but the neutron 
enters the target and is captured on Gd. The neutrino interaction vertex is 
not in the target, but there is a well measured positron event followed by a 
Gd capture signal. This interaction is thus selected as a neutrino event. 
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These two kinds of events do not compensate exactly. However, the simulation 
shows that the difference is of the order of ~1 % of the total neutrino interaction 
rate (the software used for this simulation is a low energy neutron Monte-Carlo 
that was extensively used and checked for the Bugey experiments BUG96 ). This 
imperfect compensation is due to the presence of gadolinium in the target only. 
But, since this corresponding cross section is high only at cpithermal neutron ener- 
gies, the neutrons slow down identically in both media. The difference of behavior 
happens only in the last few centimeters of the neutron path, before its capture. 
This spill in/out effect would lead to an irreducible ~1 % systematic uncertainty 
in a new single detector experiment. However, it will cancel in the Double-CHOOZ 
oscillation analysis since two identical detectors will be used. Nevertheless, a sec- 
ond order spill in/out difference will remain in Double-CHOOZ since the neutrino 
direction with respect to the neutrino target boundary changes slightly between the 
two detectors. Indeed, this small effect comes from the correlation of the V e and 
the neutron directions |CHO03| . 

In conclusion, the method used to identify a neutrino interaction allows a very 
good definition of the number of target atoms. The major concern is the design, 
the construction and the monitoring of the inner acrylic cylinders. 

4.2.2 Measurement of the fiducial volume 

We have to measure the volume of the inner acrylic vessels with an uncertainty be- 
low 0.2 %. With standard commercial materials such as dosing pumps, it is hard to 
have an absolute volume determination better than 0.5 %. We thus suggest to use 
a combination of direct and indirect measurements to obtain the required precision. 

A possible solution is to use weight measurements. For this, an intermediate 
vessel close to the acrylic target is necessary (in the experimental hall). We plan 
to measure first the weight of the empty intermediate vessel, then fill the target 
vessel and re-measure its weight. The difference of the two measurements indicates 
very accurately the mass of liquid used to fill the target. Associated with a density 
measurement, this could provide the volume measurement with uncertainty below 
0.1 % (below 10 kg on the mass determination, and around 10~ 4 on the density 
measurement). 

A second method under study consists to use pH measurement. This measure- 
ment has to be done with an acid/water mixture. It seems that this method can 
reach an 0.2 % accuracy. 

Independently of the volume of liquid used to fill the vessels, both detector have 
to be kept at the same temperature. We will thus have to monitor and control it. 
A simple regulation loop in the external veto is foreseen. 

4.3 Light collection 

We consider in the following a concentric cylindrical model of the Double-CHOOZ 
detectors consisting of the target, the 7-catcher and the outer buffer. The target 
volume of the detector is filled with organic liquid scintillator (LS) loaded with 
Gadolinium (Gd) consisting of a mixture of 

• PXE as solvent with small amount of Gd (1 g/1), 

• PPO (2,5-diphenyloxazole) as first ffuor with a concentration of 6.0 g/1, 
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• bis-MSB with concentration of 0.02 g/1 as second fluor or wavelength shifter. 

The volume of the 7-catcher enclosing the target is filled with the same LS but 
without admixture of the Gd salt. The tank containing the non scintillating buffer 
is covered by a reflecting material, and about 500 PMTs are installed on this sur- 
face (later called the PMT surface). Figure [Ol displays the PMTs mounted on the 
support structure. 

We consider the reflection coefficient (k) as a free parameter; it can be changed 
within the interval from (absolutely black surface) up to 0.98 (mirror reflection by 
VM2000 film |Mot04l l3M] L Charged particles deposit energy in the LS medium, 
mostly due to their interaction with the solvent molecules. PXE excited molecules 
transfer their energy to the PPO molecules via non-radiative processes. Then, an 
energy transfer occurs between the PPO and the bis-MSB, mainly by radiative 
transitions (100 % probability). Therefore the primary (observed) fluorescence of 
LS is connected with the radiative decay of the bis-MSB excited molecules. The 
energy spectrum of the photons emitted by the shifter is shown in Figure l4~5l The 
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Figure 4.5: Emission spectrum of the bis-MSB wavelength shifter. 

radiative transport from the light emission vertex to the PMTs is described by a 
GEANT4 Monte-Carlo simulation. Borexino-like PMTs cover between f2.5 % to 
17.5 % of the surface of the supporting cylinder. The quantum efficiency of the 
photocathode is shown on Figure POD The Monte-Carlo simulation developed for 
this work is based on the light propagation model described in jLPOOj and p3irks , 
and has been used for the Borexino experiment. The time decay of the emitted 
bis-MSB photons is described phcnomenologically by the sum of few exponentials 
having time constants ~5 ns. The light yield of the LS is taken to be 8,000 photons 
per MeV 2 , both for the target and for the 7-catcher of the detector. The photons 
emitted by the bis-MSB propagate through the target volume, and interact with 
PXE, PPO, bis-MSB and Gd salt molecules. Two physical processes have been 
taken into account: 

• (Rayleigh) elastic scattering, 

• absorption. 
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Figure 4.6: Quantum efficiency of the PMT photocathode. 



The light attenuation was described by an exponential function with the extinction 
coefficient n(X) = fi a (\) + [i s (\) , where /x a (A) is the absorption coefficient, /x s (A) the 
scattering coefficient and A the wavelength of the light. The mean free path of the 
photon is equal to A(A) = l/(log(m x /i(A))), where m is the molar concentration 
of the relevant scintillation component. The cross sections for these interactions 
have been extracted from the experimental data, obtained by usual spectroscopy 
methods. An example of A(A) variation for bis-MSB is presented in Figurc l4"?l Two 
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Figure 4.7: Absorption spectrum of the bis-MSB wavelength shifter. 



different behaviors can be seen. At wavelengths longer than 450 nm the absorbance 
drops rapidly and the measured extinction coefficients are practically equal to the 
coefficient for Rayleigh scattering, while at wavelengths shorter than 450 nm pho- 
tons absorption is the main interaction process. Elastically scattered photons have 

2 This is 2/3 of the standard unloaded pure PC scintillator. 
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an angular distribution described as 1 + cos 2 9, independent of the wavelength. The 
process of light absorption can be accompanied by an isotropic re-emission of the 
photons. The spectrum of re-emitted photons and time of the re-emission process 
were taken equal to the fluor primary spectrum and light time decay (1.3 nsec for 
bis-MSB). The re-emission probability was assumed to be equal to the absorbing 
molecule quantum efficiency taken around 0.36 for PXE, 0.8 for PPO and 0.96 for 
bis-MSB. This absorption/re-emission process can occur several times until cither 
the photon is absorbed in the scintillator volume (its energy disappears due to the 
non-radiative processes) or its wavelength falls in a region where the absorption 
probability is negligible. Photon reflection (or absorption) near the wall of the 
buffer is described by the reflection coefficient. As a result of the transport process, 
a part of the photons reaches photocathode surface of the PMTs. The spectrum 
of these photons is shown in Figure 14.81 It can be seen that the left part of the 
spectrum decreases more rapidly with respect to the emitted spectrum of bis-MSB; 
this is connected with the self-absorption of bis-MSB molecules. The results of the 
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Figure 4.8: Spectrum of the photons as they arrive at the PMT surface. 



simulation are presented as a number of photoelectrons per MeV of energy deposit, 
from point like events generated inside the target and the 7-catcher to PMT- visible 
photons that propagate to the photocathode of PMTs. If the buffer wall is black 
(reflection coefficient k=0) the number of photoelectrons was found to be around 
300 for events in the target center (for a 17.5 % coverage). This number increases 
up to 40 % if the buffer wall is reflective. The light collection time distribution is 
shown in Figure 14.91 Obviously, the reflected light increases the tail of the time 
distribution. During the first 30 ns, all photoelectrons arise from the photons that 
directly reach the PMT surface. The simulation shows a very good light collection 
homogeneity. The dependence of light collection from the event position inside the 
target was found to be within 5 % and increased up to 10 % at the position near the 
walls of the 7-catcher. The collection of the reflected light improves the homogene- 
ity and for a reflection coefficient of 0.8, the light collection is very homogeneous 
(+2-3 % at the target border, +5 % at the scintillating buffer border). 
Detector design, the number of PMTs and their positioning is now being optimized 
based on the Monte-Carlo simulation presented here. 
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Figure 4.9: Light collection for different reflectivity coefficients of the PMT support 
structure ranging from k=0 (black paint) to k=0.98 (VM2000 foil |Mot04[|gMp . 

4.4 Electronics 

4.4.1 Data recording 

The following data have to be recorded: 

• Charge and time for each PMT. 

• Pulse shape for PMT clusters, to identify recoil protons due to fast neutrons. 

Neutrino events are made of two light pulses, separated by a delay of a few fis to 
200 ^s. The single trigger rate, although lower than 10 Hz can still be reduced using 
the delayed coincidence, as explained below. This imply to store the data of the 
two pulses before the trigger decision. Furthermore, for calibration with Cf sources, 
several neutrons are detected after the fission signal, with an average multiplicity 
of 4, extending beyond 8. The dead time of the system must be kept low, stable 
and simple to control, since it contributes to systematic error. 

The front-end electronics will have to: 

• Separate the signal from the high voltage, if positive. 

• Amplify the signal by a factor ^50 to use the single photoelectron range 
(PMTs stability monitoring). 

• Add the analog signals (the total sum over the detector PMTs will be used in 
the trigger). 

• Include a discriminator per PMT (to monitor the trigger stability). 

The minimal solution for digitization is to use multihit charge ADCs, shapers and 
multihit TDCs for all channels, completed by Flash-ADCs for a few tens of PMT 
groups. The alternative would be to use Flash-ADCs for all channels, build PMT 
clusters and emulate ADCs and TDCs. A new a model of Waveform recorder with a 
smart memory management is being developed for Double-CHOOZ. It will provide 
a multihit capability and digitization with zero dead time. This is an upgrade of 
an existing model used for Borexino (the prototype will be ready in 2004). 
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4.4.2 Trigger logic 

The plan is to keep the trigger logic as simple as possible. It will be based on a 
rough energy measurement made by the analog sum of all PMT signals. A first 
level (single pulse) trigger will feature two channels: 

1. The "particle" channel: a pulse of 0.5-50 MeV, which will cause the recording 
of all channels. 

2. The "muon" channel: a pulse above 50 MeV or a signal in the veto which will 
cause the recording of time and energy information in a digital LIFO. 

The data are read out for all first level triggers and a second level trigger (final) is 
made online with the coincidence of two "particle" triggers, within 200 fis. A final 
event is composed of two singles, including information about the last muons. The 
data for each "particle" will be composed of 

• the charge and time for all PMTs, 

• the pulse shapes for ~16 PMT groups (multiplicity tunable by software). 

In addition, during data taking, some artificial light pulse patterns will be generated 
inside the target, using laser or LEDs. These artificial events will mimic the physical 
v c , in order to monitor the trigger system efficiency. Of course, each event triggered 
will carry a specific tag and serial number, for its identification in the offline data 
analysis. Table l4~D summarizes the expected rates for neutrino like triggers. The 



Table 4.1: Summary the expected trigger rates for neutrino like events at CHOOZ- 
near. Trigger rates at CHOOZ-far will be smaller. 

resulting data flow at CHOOZ-near will be around 20 kB/event, dominated by pulse 
shape data. With a trigger rate lower than 1 Hz, the amount of data remains below 
2 GB/day. 



Neutrinos 
Artificial 

Multineutron after a muon 
Cosmogenic 
Fast neutrons 
Accidental coincidence 



0.04 Hz 
0.05 Hz 
0.3 Hz 
0.001 Hz 
0.001 Hz 
0.001 Hz 
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Chapter 5 



Liquid scintillators and 
buffer liquids 

5.1 Liquid inventory 

The Double-CHOOZ detector design requires different liquids in the separate detec- 
tor volumes as shown in Figures \'6 . 41 and [3.51 The inner most volume of 12.7 m 3 , the 
z7 e -target, contains a proton rich liquid scintillator mixture loaded with gadolinium 
(Gd-LS) at a concentration of approximately 1 g/liter. The adjacent volume, the 
7-catcher, has a volume of 28 m 3 and is filled with an unloaded liquid scintillator. 
The photomultipliers are immersed in a non-scintillating buffer in order to shield 
the active volume from the gamma rays emitted by them. The volume of the buffer 
liquid is approximately 100 m 3 . Last, an instrumented volume of approximately 
110 m 3 encloses the whole setup serving as a shield against external radiation and 
as a muon veto system. Table l5~Tl summarizes the liquid inventory of a single detec- 
tor system. The selection of the organic liquids are guided by physical and technical 



Labeling 


Volume [m 3 ] 


Type 


Fe-target 


12.7 


Gd loaded LS (0.1 %) 


7-catcher 


28.1 


unloaded LS 


Buffer 


100 


non-scintillating organic liquid 


Veto 


110 


low-scintillating organic liquid 



Table 5.1: Overview of liquid inventory for a single detector. Alternatively we 
consider as well the use of a water Cherenkov detector for the veto. 

requirements, as well as by safety considerations. In particular, the solvent mixtures 
or their components have a high flash point (e.g. phenyl- xylylethane (PXE): flash 
point (fp) 145 °C, dodecane: fp 74 °C, mineral oil: fp 110 °C). The j7 c -target and 7- 
catcher have as solvent a mixture of 80 % dodecane and 20 % PXE, or alternatively 
trimethyl-benzene (PC). Mineral oil is under study as an alternative to dodecane. A 
similar solvent mixture matching the density of the 7-catcher and 77 c -target, will be 
used as the buffer liquid, however with the addition of a scintillation light quencher 
(e.g. DMP). Alternatively, pure water is under investigation provided the buoyancy 
forces can be contained, or a density matched water-alcohol mixture. The veto 
volume contains low-scintillating organic liquid and will be equipped with PMTs. 
Alternatively, we also consider to fill the veto with water and to operate it as a 
water Cherenkov detector. 
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5.2 Status of available scintillators 

Metal loading of liquid scintillators have been comprehensively studied in the frame- 
work of the LENS (Low Energy Solar Neutrino Spectroscopy) R&D phase LEN99 . 
The key groups involved in this research, MPIK and LNGS/INR, are contributing 
their expertise to the Double-CHOOZ project. The challenge of the LENS project 
was to produce stable liquid scintillators loaded with ytterbium as well as indium 
at 5-10 % in weight while simultaneously achieving attenuation lengths of several 
meters and high light yields. Novel scintillator formulations |MPI03bI |MPI04a 
ICat04al ICat04b) have been developed successfully. The scintillators have surpassed 
longterm tests on the scale of up to several years. Several prototype detectors 
filled with different scintillator samples are continuously measured in the LENS 
low-background facility at Gran Sasso since October 2003 to study the stability of 
the scintillator as well as backgrounds. No change in light yield nor in attenuation 
length has been observed and backgrounds are extremely low. 

Research with gadolinium loaded scintillator at MPIK and LNGS/INR indicates 
that suitable gadolinium loaded scintillators can be produced using the chemistry 
of beta-diketone complexes as well as using a single carboxylic acid stabilized by 
careful control of pH. Furthermore, research is being carried out to achieve stability 
with respect to interaction with detector container materials, through the adjust- 
ment of inert solvent components of the scintillator while simultaneously retaining 
high scintillation yields. 

Beta-diketonate (BDK) Gd-LS: 

The studies of the synthesis and properties of bcta-diketonates of rare earths and 
their relevant chemistry, especially stability at high temperatures, is illustrated in 
Har92, Har85 . First results of Gd-betadiketonate loaded liquid scintillators have 
been reported in |MPI03a| . Figure IB~2l displays the scintillation yield of the un- 
loaded PXE BOR04 based scintillator as a function of dodecane concentration. 
A scintillation yield of 78 % with respect to pure PXE is observed at a volume 
fraction of 80 % dodecane. Figure l5~Tl shows the light yield of a scintillator system 
with a solvent mixture of 80 % dodecane and 20 % PXE with varying PPO fluor 
concentration. The observed light yield corresponds to 80 % of the unloaded scin- 
tillator mixture, or to 60 % of a pure PC based scintillator. Attenuation length of 
the Gd-betadiketonate is being studied and values greater than 10 m at 430 nm 
have been observed after optimizing the synthesis steps. Figure IB~3l compares the 
spectral attenuation length of commercial 0.1 % Gd-acetylacetone (Gd-acac) with 
that synthesized by us. A secondary fluor(bis-MSB, emission spectrum peaked be- 
tween 420 to 450 nm) at 20-50 mg/1 is used to match the emission to the absorption 
spectrum (wavelength shifter) 

Carboxylate (CBX) Gd-LS: 

The chemical preparation of Gd loaded carboxylic acid based scintillators (single 
acid, pH controlled) has been established and demonstrated to be sound in our 
laboratories. These results have been submitted for publication, are in preparation 
for submission and are presented in publications |Cat04al ICat04bl IDan031 lMPI04b 
MPI03c . Pro gress has been swift towards the definition of scintillator specifics and 
quantitative performance. The main aspects are summarized below. 

A variety of Gd carboxylate scintillators have been produced, using methyl- 
valeric (Cq), ethyl-hexanoic (Cs) as well as trimethyl-hexanoic (Cg) acids. The 
possible solvents are trimethyl-benzene (PC) or PXE, mixed with either dodecane or 
mineral oil. The Gd scintillator can be synthesized by adding a crystalline material 
or by direct extraction into the liquid. Proper control of pH during the synthesis is 
important. 
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Figure 5.1: Scintillation light yield of 80 % dodecane 20 % PXE 0.1 % Gd beta- 
diketonate LS with varying PPO concentration relative to the unloaded 80 % do- 
decane 20 % PXE mixtures with PPO at 6 g/l. 



The solubility of two candidate Gd-carboxylate compounds namely Gd-2MeVA 
and Gd-EtHex, have been measured in a 65 % PC and 35 % Dodecane solvent 
mixture and found to be respectively 16.0 and 3.2 g/l. Light yields of 60 % with 
respect to pure PC and attenuation length of 15 m have been achieved with Gd con- 
centrations of 4 g/l and BPO (the primary fluor) concentration of 4 g/l in the same 
solvent mixture. A Cg CBX version in 50 % PC and 50 % dodecane and 1 g/l Gd 
gave 87 % of light with respect to the unloaded mixture. Good optical properties 
have been achieved. 

The first stability tests at elevated temperatures have been carried out suc- 
cessfully with the carboxylate systems. Sample mixtures of PC, mineral oil and 
Gd salt were heated to 40 °C during 18 days and mixtures with dodecane instead 
of mineral oil to 50 °C during 7 days. Figure 15/11 shows the absorption spectra of 
the PC/dodecane based Gd-carboxylate LS before and after the temperature test. 
Both the light yield and the attenuation length are stable under the test conditions. 



5.3 Scintillator definition phase 

Both the beta-diketonate and the carboxylate based Gd-LS show excellent perfor- 
mances and are viable candidate liquid scintillators for the !7 e -target. The research 
on these LS shifts now from the R&D phase to the definition phase and to qualifi- 
cation test of their use in Doublc-CHOOZ. Both Gd-LS types have to undergo long 
term tests to verify no changes in the optical performance in contact with detector 
materials. Backgrounds from radioactive trace contaminations will be studied in 
the Lens Low Background Facility (LLBF) at Gran Sasso |Mot04] . Work specific 
to the different scintillator formulation are listed below. 
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Figure 5.2: Scintillation light yield of PXE/dodecane mixture with varying dodecane 
concentration. The PPO concentration is kept constant at 6 g/1. 



BDK Gd-LS: 

The nominal BDK GD-LS candidate is based on a mixture of PXE (20 %), do- 
decane (80 %), PPO (6 g/1) and bis-MSB (50 mg/1) with a Gd loading of 0.1 % 
by weight. Future laboratory work will concentrate on further optimization of the 
chemical synthesis with special focus on questions related to the solubility and pu- 
rity of Gd-acac. The solubility has an impact on the engineering of the Gd-LS 
production scheme. Moreover, the optimization of energy transfer properties will 
be studied. A further increase in light yield by fluor optimization appears possible. 
Mineral oil (MO) will be studied in more detail as an alternative to dodecane since 
the density range of MO provides the possibility to adjust buoyancy forces applied 
to the scintillator containment vessel. A PXE (20 %) / MO (80 %) based scintillator 
can be designed matching a density in the range from 0.8 to 0.9 g/1 compared to 
0.80 g/1 for the PXE (20 %) / dodecane (80 %) mixture. 

CBX Gd-LS: 

Work on the CBX Gd-LS formulation will concentrate on the selection of the car- 
boxylic acid to use in the synthesis and on determining the chemical parameters 
relevant for the chemical stability of the solution. Possible surface induced chemical 
reactions will be investigated. Optimization of light yield and attenuation length are 
being further pursued by optimizing the synthesis as well as the solvent and fluor 
composition. The same delineations concerning solvents and densities described 
previously also apply here. 

From the results of the laboratory research, we now have two working Gd-LS 
formulations and we expect that both the BDK and CBX systems will comply with 
the design goals of Doublc-CHOOZ. The designation of the default and backup LS 
formulation will be one of the milestones during the definition phase. A further 
outcome of this phase is the detail engineering of the Gd-LS production scheme. 
This will be a critical input for the finalization of the scintillator fluid systems 
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Figure 5.3: Spectral attenuation length of Gd-acac from an optimized synthesis 
compared with a commercial purchased product. Attenuation length of approxi- 
mately 12 m is achieved at wavelength of 430 nm, corresponding to the emission 
peak of the secondary shifter. 



discussed in the next section. The final selection of the buffer and veto liquids will 
be done contingent upon the mechanical design of the containment vessels and the 
definition of the Gd-LS formulation. 

5.4 Scintillator fluid systems 

The scintillator fluid systems (SFSs) include the off-site SFS for production, pu- 
rification and storage of the Gd-LS, as well as the 7-catcher LS. A possible location 
for the off-site SFS is MPIK. The on-site SFS will be on the reactor area, close to 
the experimental location. 

The SFSs scheme envisions the production and storage of the complete Gd-LS 
for both the near and far detector, in order to assure identical proton per volume con- 
centrations. The off-site SFS will include ISO-containers for storage and subsequent 
transport to the experimental site. Moreover, it will include a purification column, 
a nitrogen purging unit, a mixing chamber, nitrogen blankets and auxiliary systems. 
A similar system, known as Module-0 Har99 , has been constructed by groups in 
this LOI associated with Borexino. Since the specifications for Module-0 are more 
demanding than required for Double-CHOOZ, no problems are anticipated. 

The on-site SFSs will consist of an area above ground close to the detector 
sites for the transport tanks which will be connected to the detector by a tubing 
system. The purpose of the on-site SFS is to transfer the different liquids from 
their transport container into the detector volumes in a safe and clean way. The 
different detector volumes will be filled simultaneously and kept at equal hydrostatic 
pressures to guarantee the integrity of the detector vessels; this will require several 
parallel lines. Details of the SFSs will be worked out during the definition phase. 
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Stability of Gd-2MeVA LS after temperature test (T-50°Cfor7 days) 
Solvent: (PC 20% tDodecane 80% ,|Gd|-4g,1, Fluors included) 
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Figure 5.4: Absorption spectra of carboxylate Gel LS prior and after temperature 
test. The sample was kept at 50 °C for 7 days. The scintillator composition consists 
of PC (20 %), dodecane (80 %), [Gd]=4 g/1 and fluors. 



Chapter 6 

Calibration 



The main goal of the calibration effort is to reach maximum sensitivity to neutrino 
oscillations by comparing the positron energy spectra measured by the CHOOZ-far 
and CHOOZ-near detectors. This is necessary for reaching the desired sensitivity to 
neutrino oscillations in Double-CHOOZ. Calculations show that a relative difference 
both in geometry (construction) and in response of detectors slightly distorts the 
ratio of the spectra in both detectors. Therefore, appropriate corrections and er- 
rors obtained on the basis of absolute and relative calibration measurements should 
be administered to the data. This should be the result of detailed Monte-Carlo 
simulations (see Chapter backed up by an extensive program of source calibra- 
tions. The calibration sources (See Table l6~T)l must be deployed regularly through- 
out the detector active volume to simulate and monitor the detector response to 
positrons, neutron captures, gammas and the backgrounds in the Double-CHOOZ 
experiment. This requires a dedicated mechanical system in order to introduce 
calibration sources into the different regions of the detector. There are a number 



Technique 


Calibrations 


Optical Fibers, Diffusive Laser ball 


Timing and Charge Slopes and Pedestals, 




attenuation length of detector components 


Neutron Sources: Am-Bc, 252 Cf 


Neutron response, relative and 




absolute efficiency, capture time 


Positron Sources: 22 Na, e8 Ge 


e + response, energy scale, trigger thresh. 


Gamma Sources: 


Energy linearity, stability, resolution, 




spatial and temporal variations. 


137 Cs 


/?", 0.662 MeV 


22 Na 


/3+, 1.275 MeV + annih 


54 Mn 


EC, 0.835 MeV 


65 Zn 


1.35 MeV 


60 Co 


EC, 1.173, 1.33 MeV 


68 Ge 


EC, (3+ 1.899 MeV + annih 


88y 


EC, 0.898, 1.836 MeV 


H neutron capture 


2.223 MeV 


241 Am- 9 Be 


(a,n) 4.44 MeV ( 12 C) 


Gd neutron capture 


Spectrum in 8 MeV window 


22S Th 


2.615 MeV 


40 K 


EC, P+,0-, 11 % 1.46 MeV 



Table 6.1: Table showing the different techniques that are available to calibrate the 
Double-CHOOZ experiment. 
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of specific tasks for a successful calibration of the detectors. These include optical 
calibrations (single photoelectron (PE) response, multiple PE response, detector 
component optical constants), electronic calibrations (trigger threshold, timing and 
charge slopes and pedestals, dead time), energy (energy scale and resolution), and 
neutron and positron detection efficiency and response. In addition, detector cal- 
ibrations must test the Monte-Carlo and analysis code to verify the accuracy of 
the simulations, throughout the detector (spatially), and during the lifetime of the 
experiment. 

6.1 Optical and electronic calibrations 

The optical calibrations are based on the experience with CHOOZ and the CTF- 
Borexino experiments. In CTF-Borexino the optical calibration consists of a UV 
pulsed-laser (jitter less than 1 ns) coupled to an optical fiber illuminating separately 
each PMT. This allows the single PE response to be measured since the amplitude 
of the pulse is tuned to approximately a single PE. This technique allows the gain, 
timing slope, charge slope and pedestals to be determined relative to individual 
PMTs and to the triggers. In addition to the optical fiber calibration, the light 
attenuation in the liquid scintillator is monitored using a diffusive laser ball source, 
as has been successfully used by the SNO experiment |SNO02j . This source illumi- 
nates all the PMTs isotropically and allows the attenuation length of the detector 
components and the PMT angular response to be measured as a function of photon 
wavelength. Finally, to ensure that we are able to veto muons with high efficiency, 
we must also calibrate the PMTs mounted on the stainless steel tank. This is done 
by also connecting optical fibers to these PMTs. The attenuation length of the 
water (or oil) shielding is measured by deploying the laser ball in this region. 

6.2 Energy calibration 

The specific signature for the detection of an electron antineutrino through inverse 
beta decay is the detection of prompt gammas from the annihilation of the positron 
and the delayed capture of a neutron several tens of [is later. While direct calibra- 
tion with an antineutrino source is impossible, it is possible to simulate each of the 
components of the antineutrino signal, such as the prompt positron and delayed 
neutron by deploying positron, neutron, and gamma sources. 

The standard calibration system will include a permanent vertical tube, entering 
the detector until the center of the inner acrylic target. This open tube will allow 
frequent and safe calibration with radioactive sources. 

6.2.1 Gamma ray sources 

Positron annihilates at rest and produces 2 back-to-back gammas. Thus, for a high 
detection efficiency we must be able to calibrate the detector energy response to 
gammas from I MeV to ~10 MeV corresponding to the endpoint of the fission 
product beta decays. In addition, a neutron is detected by its capture on the Gd 
additive to the liquid scintillator and produces a gamma cascade of approximately 
8 MeV. For this reason, it is necessary to also know the energy scale in the high 
energy window of 6-10 MeV to be able to identify the delayed second trigger as a 
neutron. Specifically, it will be necessary to know the gamma energy correspond- 
ing to the neutron detection threshold for both the near and far detector with a 
100 keV accuracy. This is accomplished by deploying various higher energy gamma 
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calibration sources (see Table |6~T|) and by detailed Monte-Carlo simulations in the 
energy region where there are no calibration sources. 

The overall energy scale can be determined from the position of the 0.662 MeV 
peak of the 137 Cs source, and then verified by calibration with several gamma 
sources (see TableED in different energy ranges: 54 Mn (0.835 MeV), 22 Na (1.275 MeV), 
65 Zn (1.351 MeV), 60 Co, and 228 Th (2.614 MeV). These gammas allow the en- 
ergy response to the positron annihilation photons to be determined for different 
positron energies. The capture of neutrons from an Am-Be source scintillator (to 
be discussed later) can also be used as a high energy gamma source as it produces 
prompt 4.4 MeV gammas. We will also use the natural sources from radioactive 
impurities of the detector materials ( 40 K, 208 T1 ) and neutrons produced by cosmic 
muons for energy calibration. Since these sources are present permanently, they 
are useful for monitoring the stability of the energy response. Thus, the primary 
purpose of the gamma sources are to determine and monitor the energy scale for 
both the far and near detectors as a function of position and time during which the 
experiment is conducted. 

6.2.2 Positron response 

Positron detection can be simulated experimentally by means of the 22 Na source. 
A 22 Na source emits a 1.275 MeV primary gamma accompanied by a low energy 
positron which annihilates inside the source container. The primary and annihila- 
tion gammas from the source mimic the positron annihilation resulting from an an- 
tineutrino event inside the detector. An alternative positron source is a 68 Ge source 
which produces positrons with higher energies, and therefore calibrates higher en- 
ergy positrons. 68 Ge decays by EC to 68 Ga and /3 + -decays to stable 68 Zn with an 
endpoint of 1.9 MeV. This isotope also has the advantage that it produces only low 
energy gammas in coincidence with the nuclear decay, and the (3 + has an endpoint 
of 1.889 MeV 89 % of the time. A second purpose of this source (if a source is 
constructed so that the beta is absorbed by the shielding surrounding the source) is 
to tune the trigger threshold to be sensitive to annihilation gammas and to monitor 
its stability. A 68 Ge source has been successfully used in the Palo Verde reactor 
neutrino experiment |PV97| . 

6.2.3 Neutron response 

Coincident with the production of a positron in inverse beta decay, a neutron is 
produced. The neutron then quickly thermalizes and is captured on the Gd ( 155 Gd 
or 157 Gd, with cross sections of 60,900 and 254,000 barns, respectively) loaded in the 
central target. The neutron capture is accompanied by the emission of a cascade 
of gamma-rays with the summed energies of 8.536 and 7.937 MeV, respectively. 
Thus, neutrons are selected by cutting on gammas with energies exceeding 6 MeV. 
However, a fraction of the gamma-rays can escape detection, especially events that 
occur near the boundary of the fiducial volume. Therefore, it is expected that 
the neutron detection efficiency decreases for events near the borders of the acrylic 
vessel that contains the Gd loaded liquid scintillator. Calibration of this effect must 
be quantified by deployment of neutron calibration sources throughout the detector 
and comparing the detector response to Monte-Carlo. In addition to measuring the 
neutron response, neutron calibration is also a very sensitive method for determining 
in-situ various liquid scintillator properties, such as the hydrogen and gadolinium 
concentration in the liquid scintillator. 

There are two suitable and accessible neutron sources for neutron calibration: 
the Am-Be source and 252 Cf spontaneous fission source. These sources emit neu- 
trons with different energy spectra from what is expected from inverse beta decay, 
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and thus the importance of these differences needs to be quantified. To decrease 
the background during neutron source deployment, neutrons from Am-Be should 
be tagged by the 4.4 MeV gamma emitted in coincidence with the neutron. This 
allows the neutron capture detection efficiency to be determined independent of 
knowing the precise rate of the neutron source, because every time a 4.44 MeV 
gamma is detected a neutron is released |Cro89| . The absolute neutron detection 
efficiency can also be determined with a 252 Cf source by using the known neutron 
multiplicity (known to 0.3 %). For the source placed into the center, the size of the 
Gd region is larger than the neutron capture mean free path, so that the neutron 
capture is studied independent of the presence of the acrylic vessel. In order to tag 
the neutron events, a small fission chamber is used to detect the fission products. 
Therefore, neutron source calibrations provide us with the relevant data to cali- 
brate the detector response to neutrons. In particular, neutron sources allow us to 
measure the absolute neutron efficiency, to determine and monitor the appropriate 
thresholds of neutron detection, and to measure the neutron capture time for both 
the far and near detectors. 

6.2.4 The Calibration source deployment system 

A mechanical system must be developed to introduce calibration sources throughout 
the detector active volume. The system must be easy to set up so that calibration 
can be done frequently without loss of neutrino live time. The suggestion is to use 
a system of ropes and pulleys similar to the SNO experiment (see Figure IfTTl for a 
conceptual design). However, unlike the SNO experiment we must be able to deploy 
sources throughout the active volume, rather than in a plane as is done in the SNO 
experiment. The reason for this is that because during the lifetime of the detector, 
PMT mortality might result in an anisotropy in the detector response. Moreover, 
the effect will manifest as a anisotropy relative to Chooz-Near and Chooz-Far which 
will impact on the energy resolution and scale of the two detectors. The system of 
ropes and pulleys must be designed so that the calibration sources sample a large 
fraction of the active volume and can calibrate this effect. 

6.2.5 Map of the Gd-LS target 

The starting point of the design is to introduce sources through a glove box situated 
at the center and top of the cylinder housing the Gd-LS target. In the glove box 
sources can be prepared for deployment without introducing contaminants into the 
active volume. The glove box can also be evacuated and flushed with LN2 before 
deploying sources to prevent Radon from entering the active volume. The sources 
are then suspended from a rope and lowered straight down from the glove box to 
the bottom of the cylinder using a stepper motor. In this way we can calibrate 
the variation of the detector response along the axis of symmetry of the cylinder 
(z-axis). To calibrate the detector response off of the z-axis the calibration sources 
must be physically moved away from the z-axis such that the detector response as 
a function of radius can be determined. To achieve this, the idea is to fasten two 
ropes on opposite sides of the cylinder. Then, to feed the rope through two pulleys 
(one for each side rope) attached to top of the calibration sources, then to a stepper 
motor located near the glove box. The tension on either of the ropes can then be 
independently adjusted by carefully controlling stepper motors. When the tension 
is changed of one of the side ropes the source will move away from the z-axis. This 
will allow the sources to be deployed throughout most of the area of a plane defined 
by the central axis of symmetry (z-axis) and the line connecting the places that the 
ropes are attached. Including as second set of ropes perpendicular to the first will 
allow the source to be moved not only away from the z-axis but throughout most 
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of the active volume of the target. The SNO experiment has been able to attain 
a deployment accuracy of 5 cm using this method when the source is moved in a 
plane [NIMOlj . 

6.2.6 Calibrating the gamma-catcher, buffer, and veto 

A deployment mechanism must also be devised to deploy sources outside of the 
main central target. Sources must be deployed in the gamma catcher region as well 
as the buffer and the veto. A further requirement of the components inside the 
active volume is that the system must not block the scintillation light, nor change 
the detector response, and they must be impervious to liquids. The suggestion 
is again to use a system of pulleys and side ropes to cover most of the volume. A 
possible scenario for such a system is shown in Figure lETTl The system would sample 
calibration source positions in a plane from the z-axis outward through the cylinders 
to the veto cylinder. The sources will be accessed through a glove box which will 
be movable so that it can be mounted on top of all the source interfaces. The 
calibration of the veto can be done with a rail deployment system (only 1 central 
rope but movable along the radius of the cylinder) , since here the mechanism can be 
constructed without blocking the scintillation light. The right panel of Figure l6~Tl 
shows a possible configuration of the ropes, specifically the figure shows how the 
top (ropes outward from the center) and sides (ropes in shape of squares) of the 
detector will be calibrated. However, calibration of the bottom of both detectors 
is more difficult since access to this region is limited. Calibration of the bottom 
portion of the detectors still needs to be investigated. 



Chapter 7 

Backgrounds 



The signature for a neutrino event is a prompt signal with a minimal energy of 
about 1 MeV and a delayed 8 MeV signal after neutron capture in gadolinium. 
This may be mimicked by background events which can be divided into two classes: 
accidental and correlated events. The former are realized when a neutron like event 
by chance falls into the time window (typically few 100 /zs) after an event in the 
scintillator with and energy of more than one MeV. The latter is formed by neu- 
trons which slow down by scattering in the scintillator, deposit > 1 MeV visible 
energy and are captured in the Gd region. In this chapter we first discuss possible 
sources and fluxes for background events and later estimate their rates. With these 
numbers we find criteria for the necessary overburden of the near detector and we 
will extract purity limits for detector components. 



7.1 Beta and gamma background 
7.1.1 Intrinsic beta and gamma background 

In this section the intrinsic background due to beta and gamma events above 
~1 MeV is discussed. It can be produced in the scintillator or in the acrylic vessels 
which contain the liquid. The contribution from the Uranium and Thorium chains 
is reduced to a few elements, as all alpha events show quenching with visible ener- 
gies well below 1 MeV. Furthermore the short delayed Bi-Po coincidences in both 
chains can be detected event by event, and hence rejected. In the end, only the 
decays of 234 Pa (beta decay, Q = 2.2 MeV), 228 Ac (beta decay, Q = 2.13 MeV) and 
208 T1 (beta decay, Q — 4.99 MeV) have to be considered. Assuming radioactive 
equilibrium the beta/gamma background rate due to both chains can be estimated 
by h ~ Mu • 6 • 10 3 s^ 1 + M Th • 4 • 10 3 s _1 , where the total mass of U and Th 
is given in gram. Taking into account the total scintillator mass of the neutrino 
target plus the 7-catcher, this rate can be expressed by b\ ~ 3 s~ 1 (cu,Th/10 -11 ) , 
where cu.Th is the mass concentration of Uranium and Thorium in the liquid. The 
contribution from 40 K can be expressed by 62 — 1 s _1 (ck/10~ 9 ), where ck is the 
mass concentration of natural K in the liquid. 

The background contribution due to U, Th and K in the acrylic vessels can 
be written as 63 ~ 2 s _1 (aK/10 -7 ) + 5 s _1 (au,Th/10 -9 ), where ax and au.Th de- 
scribe the mass concentrations of K, U and Th in the acrylic. In total, the intrinsic 
beta/gamma rate is the sum b = b\ + 62 + &3- In the CTF of the Borexino exper- 
iment at Gran Sasso, concentration values of cu,Th < 10~ 15 and ck < 10~ 12 have 
been measured for two liquid scintillators (PC and PXE) with volumes of about 
4 m 3 . Upper limits on radioactive trace elements in acrylic have been reported 
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to be au,Th < 3 ■ 10~ 12 by the SNO collaboration SNO02 . Gamma spectroscopy 
measurements show upper limits of ok < 1 • 10~ 9 . This shows that in principle the 
beta/gamma rate in the detector due to intrinsic radioactive elements can be kept 
at levels well below 1 s _1 . The aimed concentration values for this goal are given 
in Table O 



Element 


allowed concentration (g/g) 




for b < 1 s _1 


Uranium, Thorium in scintillator 


~ 10"" 


Potassium in scintillator 


~ io- 10 


Uranium, Thorium in acrylic vessels 


~ 1(T 10 


Potassium in acrylic vessels 


~ IO" 8 



Table 7.1: Upper limits on U, Th and K concentrations in the liquid scintillator 
and acrylic vessels to achieve a beta/gamma rate below 1 s _1 



7.1.2 External gamma background 

According to the experience gained in the CTF of Borexino the dominant contribu- 
tion to the external gamma background is expected to come from the photomulti- 
pliers (PMTs) and structure material. Again contributions from U, Th and K have 
to be considered. However, because of the shielding of the buffer region only the 
2.6 MeV gamma emission from 208 T1 has to be taken into account. The activity of 
one PMT in the CTF (structure material included) is known to be ~ 0.4 s _1 . The 
shielding factor S due to the buffer liquid can be calculated to be S ~ 10~ 2 . Hence, 
the resulting gamma background in the neutrino target plus the 7-catcher can be 
written as b ext — 2 s -1 (-/Vpmt/500), where Npmt is the number of PMTs. 

7.2 Neutron background 

7.2.1 Intrinsic background sources 

Neutrons inside the target may be produced by spontaneous fission of heavy ele- 
ments and by (a,n)-reactions. For the rate of both contributions the concentrations 
of U and Th in the liquid are the relevant parameters. The neutron rate in the 
target region can be written as n int — 0.4 s _1 (cu,Th/10 -6 ). Hence, for the aimed 
concentration values as described above the intrinsic contribution to the neutron 
background is negligible. 

7.2.2 External background sources 

Several sources contribute to the external neutron background. We first discuss 
external cosmic muons which produce neutrons in the target region via spallation 
and muon capture. Those muons intersect the detector and should be identified by 
the veto. However, some neutrons may be captured after the veto time window. 
Therefore we estimate the rate of neutrons, which are generated by spallation pro- 
cesses of through going muons and by stopped negative muons which are captured 
by nuclei. 

The first contribution is estimated by calculating the muon flux for different 
shielding values and taking into account a E - 75 dependence for the cross section of 
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neutron production, where E is the depth dependent mean energy of the total muon 
flux. The absolute neutron flux is finally obtained by considering measured values 
in several experiments (LVD |LVD99| . MACRO )MA(J98| . CTF |0TF98| 1 in the 
Gran Sasso underground laboratory and extrapolating these results by comparing 
muon fluxes and mean energies for the different shielding factors. Table [7~2"1 gives 
the expected neutron rate depending on the shielding. 



Overburden 


Muon rate 


Mean muon energy 


Neutrons 


(m.w.e.) 


(s- 1 ) 


(GeV) 




40 


1.1 • 10 ;i 


14 


2 


60 


5.7- 10 2 


19 


1.4 


80 


3.5 • 10 2 


23 


1 


100 


2.4- 10 2 


2G 


0.7 


300 


2.4- 10 1 


G3 


0.15 



Table 7.2: Estimated neutron rate in the active detector region due to through 
going cosmic muons. 

Negative muons which are stopped in the target region can be captured by nuclei 
where a neutron is released afterwards. The rate can be estimated quite accurately 
by calculating the rate of stopped muons as a function of the depth of shielding and 
taking into account the ratio between the /x-life time and /^-capture times. As the 
capture time in Carbon is known to be around 25 /is («1 ms in H) only about 10 % 
of captured muons may create a neutron. Since the concentration in Gd is so low, 
its effect can be neglected here. The estimated results are shown in Table 1731 The 
neutron generation due to through going muons dominates. 



Overburden 


Muon stopping rate 


Neutrons 


(m.w.e.) 




(s- 1 ) 


40 


5 • 10" 1 


0.7 


GO 


3 • IO" 1 


0.4 


80 


1.2 • 10" 1 


0.2 


100 


6 • 10" 2 


0.08 


300 


2.5 • 10" 3 


0.003 



Table 7.3: Estimated neutron rate in the target region due to stopped negative 
muons. 



7.2.3 Beta-neutron cascades 

Muon spallation on 12 C nuclei in the organic liquid scintillator may generate 8 He, 
9 Li, and 11 Li which may undergo beta decay with a neutron emission. In that case 
those background events show the same signature as a neutrino event. For shallow 
shielding depths the muon flux is too high to allow tagging by the muon veto, as 
the lifetimes of these isotopes are between 0.1 s and 1 s. The cross sections for the 
production of 8 He, 9 Li have been measured by a group of TUM at the SPS at CERN 
with muon energies of 190 GeV (NA54 experiment HagOO ). In this experiment only 
the combined production 8 He + 9 Li where obtained without ability to separate each 
isotope. An estimate for the background rates for shallow depth experiments like 
CHOOZ can be obtained from results of the KamLAND experiment by calculating 
the muon flux for energies above ~500 GeV |Hor03| . With this assumption an event 
rate of about 0.4 per day in the target region can be estimated for a 300 m.w.e. 



8He 119ms 




&Be 



2 4He 

a : 0.09 MeV 

Figure 7.1: Relevant branching ratios for the decay of the 8 He isotope, normalized 
to 100 %. Half-lives are quoted, as well as the end-point of the (3 decays. Neutrons 
emitted in these decays are typically around 1 MeV. The double cascade decay to 
the 8 Be offer a possibility to measure , in situ, the production rate. 



shielding. A more conservative estimate is obtained assuming a E 75 scaling as we 
did in calculating the neutron flux. Then the rate should be around 2 events per 
day. In Table 17^1 all radioactive 12 C-spallation products including the beta-neutron 
cascades are shown with the estimated event rates in both detectors. 

The Q-values of the beta-neutron cascade decays is 8.6 MeV, 11.9 MeV, 20.1 MeV 
for 8 He, 9 Li, and 11 Li, respectively. In the experiment the 8 He production rate might 
be measured if we set a dedicated trigger after a muon event in the target region 
looking for the double cascade of energetic betas ( 8 He — > 8 Li — > 8 Be) occuring in 
50 % of all 8 He decays (see Figures 17. II and l7T2*|l . Nothing similar exist in the case 
of the 9 Li, but the beta endpoint is here above the endpoint of positron induced by 
reactor antineutrinos. Nevertheless, from the NA54 experiment |Hag00| results the 
total cross section of 8 Hc + 9 Li is known, and if the 8 He is evaluated separately, 
some redundancy on the total /3-neutron cascade will be available. Figure mi shows 
the relevant branching ratios of the 8 He isotope, normalized to 100 %. The neutrons 
emitted in these decays are typically around 1 MeV. Figure f7~2l shows the relevant 
branching ratios of the 9 Li isotope, normalized to 100 %. 

7.2.4 External neutrons and correlated events 

Very fast neutrons, generated by cosmic muons outside the detector, may penetrate 
into the target region. As the neutrons are slowed down through scattering, recoil 
protons may give rise to a visible signal in the detector. This is followed by a de- 
layed neutron capture event. Therefore, this type of background signal gives the 
right time correlation and can mimic a neutrino event. Pulse shape discrimination 




Figure 7.2: Relevant branching ratios for the decay of the 9 Li isotope, normalized 
to 100 %. Half-lives are quoted, as well as the end-point of the j3 decays. Neutrons 
emitted in these decays are typically around 1 MeV. In case of (3 decay to 8 Be, the 
latter transform immediately to two low energy a particles. 



in order to distinguish between (3 events and recoil protons is in principle possible, 
but should not be applied in the analysis as additional statistical and systematic 
uncertainties should be avoided in the experiment. As the muon is not seen by the 
veto, those correlated events may be dangerous for the experiment. 

Therefore a Monte-Carlo program has been written to estimate the correlated 
background rate for a shielding depth of 100 m.w.e. and flat topology. In order 
to test the code the correlated background for the old Chooz experiment (different 
detector dimensions, 300 m.w.e. shielding) has been calculated with the same pro- 
gram. The most probable background rate was determined to be 0.8 counts per day. 
A background rate higher than 1.6 events per day is excluded by 90 % C.L. This 
has to be compared with the measured rate of 1.1 events per day. We conclude that 
the Monte-Carlo program reproduces the real correlated background value within 
roughly a factor 2. 

For Double-Chooz we calculated the correlated background rate for 100 m.w.e. 
shielding and estimated the rates for other shielding values by taking into account 
the different muon fluxes and assuming a E 0,75 scaling law for the probability to 
produce neutrons. The neutron capture rate in the Gd loaded scintillator for an 
overburden of 100 m.w.e. is about 300/h. However, only 0.5 % from those neutrons 
create a signal in the scintillator within the neutrino window (i.e. between 1 MeV 
and 8 MeV), because most deposit in total much more energy during the multiple 
scattering processes until they are slowed down to thermal energies. The quenching 
factors for recoil protons and carbon nuclei has been taken into account. In addition 
around 75 % from those events generate a signal in the muon veto above 4 MeV 
(visible [3 equivalent energy). In total the correlated background rate is estimated 
to be about 3.0 counts per day for 100 m.w.e. shielding. In Table l7~5l the estimated 
correlated background rates are shown for different shielding depths. 

The correlated background rates can be compared with accidental rates, where 
by chance a neutron signal falls into the time window opened by a /3 + -like event. 
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9 Li 
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3.6 


1.7 ±0.3 


0.36 


8 Li 


31 ±12 


6.6 


3.3 ± 1.2 
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8 He 




8 He & 9 Li measured together 




6 He 


126 ±12 


26.8 


13.2 ± 1.3 


2.8 


n c 


7100 ± 455 


1510 


749 ± 48 


159.3 


io c 


904 ± 114 


192 


95 ± 12 


20.2 


9 C 


38 ± 12 


8.1 


4.0 ± 1.2 


0.85 


8 B 


60 ±11 


12.7 


5.9 ± 1.2 


1.25 


7 Be 


1800 ± 180 


382.9 


190 ±19 


40.4 



Table 7.4: Radioactive isotopes produced by muons and their secondary shower 
particles in liquid scintillator targets at the CHOOZ near and far detectors. The 
rates (events/d) are given for a target of 4.4 x 10 29 12 C (For a mixture of 
80 % Dodecane and 20 % PXE, 12.7 m 3 ) at a depth of 60 m.w.e. for the near 
detector and 300 m.w.e. for the far detector. Because of the positron annihilation 
the visible energy in j3 + decays is shifted by 1.022 MeV. 9 Li and 8 He could not be 
evaluated separately. Columns 3 and 5 correspond to an estimate of the number of 
events assuming that the isotopes are produced only by high energy muon showers 
E > 500 GeV |Hor03j . A neutrino signal rate of 85 events per day is expected at 
CHOOZ-far, without oscillation effect (for a power plant running at nominal power, 
both dead time and detector efficiency are not taken into account here). 



The background contribution due to accidental delayed coincidences can be deter- 
mined in situ by measuring the single counting rates of neutron-like and /3 + -like 
events. Therefore the accidentals are not so dangerous as correlated background 
events. Taking for granted we reach reasonably low concentrations of radioactive 
elements in the detector materials, especially in the scintillator itself (see discussion 
above), the beta-gamma rate above 1 MeV can be expected to be about a few counts 
per second. If the time window for the delayed coincidence is ^200 /is (this should 
allow a highly efficient neutron detection in Gd loaded scintillators), and the veto 
efficiency is at 98 % the accidental background rates can be estimated as depicted 
in Table 17.61 The rate of neutrons which cannot be correlated to muons ( "effective 
neutron rate") is calculated by n c g — n tot ■ (1 — e), where n to t is the total neutron 
rate (sum of the numbers given in Table [7"21 and IT3|l and e is the veto efficiency. If 
the veto efficiency is 98 % or better, the accidental background for the far detector 
is far below one event per day (see following Table ITfi)) . 



7.2.5 Conclusion 

We conclude that correlated events are the most severe background source for the 
experiment. Two processes mainly contribute: /3-neutron cascades and very fast 
external neutrons. Both types of events are coming from spallation processes of 
high energy muons. In total the background rates for the near detector will be 
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40 


829 
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400 
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57 


0.5 



Table 7.5: Estimated neutron rate in the target region and the correlated back- 
ground rate due to fast neutrons generated outside the detector by cosmic muons. 



Overburden 


Effective neutron rate 


Accidental background rate 


(m.w.e.) 


(0 


(O 


40 


97 


2.4 


60 


65 


1.6 


80 


43 


1.0 


100 


28 


0.7 


300 


6 


0.15 



Table 7.6: Example of estimated accidental event rates for different shielding 
depths. The rates scale with the total beta-gamma rate above 1 MeV (here 
btot = b ex t + b w 2.5 s _1 ), the time window (here r = 200 /is) and the effective 
neutron background rate (here a muon veto efficiency of 98 % was assumed). 



between 9/d and 23/d if a shielding of 60 m.w.e. is choosen. For the far detector a 
total background rate between 1/d and 2/d can be estimated. 



Chapter 8 



Experimental Errors 



8.1 From CHOOZ to Double-CHOOZ 

In the first CHOOZ experiment, the total systematic error amounted to 2.7 %. The 
goal of Double-CHOOZ is to reduce the overall systematic uncertainty to 0.6 %. 
A summary of the CHOOZ systematic errors is given in Table 18.11 CHO03 . The 
right column presents the new experiment goals. Lines 1,4, and 5 correspond to sys- 
tematic uncertainties related to the reactor flux and the cross section of neutrinos 
on the target protons. These errors become negligible if one uses two antineutrino 
detectors located at different baselines. In order to improve the systematic uncer- 
tainties related to the detector and to the V e selection cuts, the Double-CHOOZ 
experiment will take advantage of the latest technical developments achieved by the 
recent sci ntillator detector CHOOZ |CHO03j . CTF |CTF98| . KamLAND |KAM02j . 
Borexino |Sch99j . and the LENS R&D phase |LEN99j . 



CHOOZ Double-CHOOZ 



Reactor cross section 


1.9 


% 




Number of protons 


0.8 


% 


0.2 % 


Detector efficiency 


1.5 


% 


0.5 % 


Reactor power 


0.7 


% 




Energy per fission 


0.6 


% 





Table 8.1: Overview of the systematic errors of the CHOOZ and Double-CHOOZ 
experiment. 



8.2 Relative normalization of the two detectors 

The goal of Double-CHOOZ is to use two v a detectors in order to cancel or decrease 
significantly the systematic uncertainties that limit the #13 neutrino mixing angle 
measurement. However, beside those uncertainties, the relative normalization be- 
tween the two detectors is the most important source of error and must be carefully 
controlled. This section covers the uncertainties related to the V c interaction and 
selection in the analysis, as well as the electronics and data acquisition dead times. 
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8.3 Detector systematic uncertainties 



8.3.1 Solid angle 

The distance from the CHOOZ detector to the cores of the nuclear plant have been 
measured to within ±10 cm by the CHOOZ experiment. This translates into a 
systematic error of 0.15 % in Double-CHOOZ, because the effect becomes relatively 
more important for the near detector located 100-200 meters away from the reactor. 
Specific studies are currently ongoing to guarantee this 10 cm error. Furthermore, 
the "barycenter" of the neutrino emission in the reactor core must be monitored 
with the same precision. In a previous experiment at Bugey BUG96 , a 5 cm change 
of this barycenter was measured and monitored, using the instrumentation of the 
nuclear power plant [Har92 . Our goal is to confirm that this error could be kept 
below 0.2 %. 

8.3.2 Number of free protons in the target 

8.3.2.1 Volume measurement 

In the first CHOOZ experiment, the volume measurement was done with an absolute 
precision of 0.3 % |CHO03| . The goal is to reduce this uncertainty by a factor of 
two, but only on the relative volume measurement between the two inner acrylic 
vessels (the other volumes do not constitute the v c target). An R&D program 
has already started in order to find the optimal solution for the relative volume 
determination (See Section I4.2.2|) . Among some ideas under study, we plan to 
use the same mobile tank to fill both targets; a pH-based measurement is being 
studied as well. A more accurate measurement could be performed by combining a 
traditional flux measurement with a weight measurement of the quantity of liquid 
entering the acrylic vessel. Furthermore we plan to build both inner acrylic targets 
at the manufacturer and to move each of them as a single unit into the detector site. 
A very precise calibration of both inner vessels is thus foreseen at the manufacturer 
(filling tests). 

8.3.2.2 Density 

The uncertainty of the density of the scintillator is ^0.1 %. The target liquid will 
be prepared in a large single batch, so that they can be used for the two detector 
fillings. The same systematic effect will then occur in both detectors and will not 
contribute to the overall systematic error (this effect will be included automatically 
in the absolute normalization error, see Chapter [[J}. However, the measurement 
and control of the temperature will be mandatory to guarantee the stability of the 
density in both targets (otherwise it would contribute to the relative uncertainty, see 
Chapter |5J). To thermalize both 77 e targets, the temperature control and circulation 
of the liquid in the external veto is foreseen. 

8.3.2.3 Number of hydrogen atoms per gramme 

This quantity is very difficult to measure, and the error is of the order of 1 %; 
however, the target liquid will be prepared in a large single batch (see above). This 
will guarantee that, even if the absolute value is not known to a high precision, 
both detectors will have the same number of hydrogen atoms per gramme. This 
uncertainty, which originates in the presence of unknown chemical compounds in 
the liquid, does not change with time. 



8.3.3 Neutron efficiency 



The thermal neutron is captured either on hydrogen or on Gadolinum (other reac- 
tions such as Carbon captures can be neglected). We outline here the systematical 
errors related to the neutron signal. 

8.3.3.1 Gadolinium concentration 

Gd concentration can be extracted from a time capture measurement done with a 
neutron source calibration (see Chapter |SJ. A very high precision can be reached 
on the neutron efficiency (0.3 %) by measuring the detected neutron multiplicity 
from a Californium source (Cf). This number is based on the precision quoted 
in |CM(J03| , but taking away the Monte-Carlo uncertainty, since we work with two- 
identical detectors. This precision is expected to be better by a factor of two in 
the Double-CHOOZ experiment because it is easier to compare two experimental 
measurements in identical detectors than to compare a theoretical spectrum with 
a measurement. We can increase our sensitivity to very small differences in the 
response from both detectors by using the same calibration source for the mea- 
surements. The Californium source calibration can be made all along the z-axis of 
the detector, and is thus snsitive to spatial effects due to the variation of Gd con- 
centration (staying far enough from the boundary of the target, and searching for 
a top/down assymetry). A difference between the time capture of both detectors 
could also be detected with a sensitivity slightly less than 0.3 %. 

8.3.3.2 Spatial effects 

We consider here the spill in/out effect, i.e the edge effect associated with neutron 
capture close to the acrylic vessel surrounding the inner target [CHO03] . and the 
angle between the neutron direction and the edge of the acrylic target that is slightly 
different between the two detectors. The ~1 % spill in/out effect oberved in the 
first CHOOZ experiment |CH()03| cancels by using a set of two identical detectors 
(same effect). Nevertheless the second effect (angle) persists, but is considered to 
be negligible. 

8.3.4 Positron efficiency 

The simulation of the Double-CHOOZ detectors confirms that a 500 keV energy 
cut induces a positron inefficiency smaller than 0.1 % (see Chapter 0J. The relative 
uncertainties between both detectors lead thus to an even smaller systematic error 
and is therefore negligible. 

8.4 Selection cuts uncertainties 

The analysis cuts are potentially important sources of systematic errors. In the first 
CHOOZ experiment, this amounted in total to 1.5 % [CHOPS] . The goal of the new 
experiment is to reduce this error by a factor of three. The CHOOZ experiment used 
7 analysis cuts to select the F e (one of them had 3 cases, see Section 8.7 of |CHO03p . 
In Double-CHOOZ we plan to reduce the number of selection cuts to 3 (one of them 
will be very loose, and may not even be used). This can be achieved because of 
reduction of the number of accidentals background events, only possible with the 
new detector design (see Chapter . To select v e events we have to identify the 
prompt positron followed by the delayed neutron (delayed in time and separated in 
space). The trigger will require two local energy depositions of more than 500 keV 
in less than 200 fis. 



8.4.1 Identifying the prompt positron signal 

Since any V e interaction deposits at least 1 MeV (slightly less due to the energy 
resolution effect) the energy cut at 500 keV does not reject any V e events. As a 
consequence, there will not be any systematic error associated with this cut (see 
Figure l%TTJ) . The only requirement is the stability of the energy selection cut, which 
is related to the energy calibration (see Chapter EJ). 
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Figure 8.1: Simulation of the positron energy spectrum (in MeV) measured with the 
Double-CHOOZ detector (10,000 events, without backgrounds). Positron energy 
is fully contained with a probability of 99.9 %, as a consequence of the 60 cm 
scintillating buffer. 



8.4.2 Identifying the neutron delayed signal 

The energy spectrum of a neutron capture has two peaks, the first peak at 2.2 MeV 
tagging the neutron capture on hydrogen, and the second peak at around 8 MeV 
tagging the neutron capture on Gd (see Figure l8~2*|l . The selection cut that identifies 
the neutron will be set at about 6 MeV, which is above the energy of neutron capture 
on hydrogen and all radioactive contamination. At this energy of 6 MeV, an error 
of ^100 keV on the selection cut changes the number of neutrons by ^0.2 %. This 
error on the relative calibration is achievable by using identical Cf calibration source 
for both detectors (see Chapter 0). 

8.4.3 Time correlation 

The neutron time capture on Gd in the CHOOZ detector is displayed in Figure 
But since the exact analytical behaviour describing the neutron capture time on 
Gd is not known, the absolute systematic error for a single detector cannot be 
significantly improved with respect to CHOOZ CHO03 . However, the uncertainty 
originating from the liquid properties disappears by comparing the near and far 
detector neutron time capture. The remaining effect deals with the control of the 
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Figure 8.2: Simulation of the neutron energy spectrum (in MeV) measured with 
the Double-CHOOZ detector (10,000 events, without backgrounds). There are two 
energy peaks for the neutron capture on hydrogen (releasing 2.2 MeV) and on 
gadolinium (releasing about 8 MeV). The Double-CHOOZ experiment will select 
all neutron events with an energy greater than 6 MeV. The resulting systematic 
uncertainty thus depends on the relative calibration between the near and far de- 
tectors. 



electronic time cuts. For completeness, a redundant system will be designed in order 
to control perfectly these selection cuts (for example time tagging in a specialized 
unit and using Flash- ADC's). 

8.4.4 Space correlation 

The distance cut systematic error (distance between prompt and delayed events) was 
published as 0.3 % in the CHOOZ experiment |CHO03| . This cut is very difficult to 
calibrate, since the rejected events are typically F e candidates badly reconstructed. 
In Double-CHOOZ, this cut will be either largely relaxed (two meters instead of 
one meter for instance) or totally suppressed, if the accidentals event rate is low 
enough, as expected from current simulations (see Chapter [JJ. 

8.4.5 Veto and dead time 

The Double-CHOOZ veto will consist of a liquid scintillator and have a thickness of 
60 cm liquid scintillator at the far site, and even larger at the near detector site. The 
veto inefficiency comes from the through going cables and the supporting structure 
material. This inefficiency was low enough in the first experiment, and should be 
acceptable for the CHOOZ-far detector. However, it must be lowered for the near 
detector because the muon flux is a factor 30 higher for a shallower overburden of 
60 m.w.e.. A constant dead time will be applied in coincidence with each through 
going muon. This has to be measured very carefully since the resulting dead time 
will be very different for the two detectors: a few percent at the far detector, and 
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Figure 8.3: Neutron delay distribution measured with the Am/Be source at the 
detector centre in the CHOOZ detector |CHO03j . The time origin is defined by the 
4.4 MeV 7-ray. 



at moreless 30 % at the near detector. A 1 % precision on the knowledge of this 
dead time is mandatory. This will require the use of several independent methods: 

• the use of a synchronous clock, to which the veto will be applied, 

• a measurement of the veto gate with a dedicated flash ADC, 

• the use of an asynchronous clock that randomly generates two particles mim- 
icking the antincutrino tag (with the time between them characteristic of the 
neutron capture on Gd). With this method, all dead times (originating from 
the veto as well as from the data acquisition system) will be measured si- 
multaneously. The acquisition of a few thousands such events per day would 
achieve the required precision, 

• the generation of sequences of veto-like test pulses (to compare the one pre- 
dicted dead time to the actually measured). 

8.4.6 Electronics and acquisition 

The trigger will be rather simple. It will use only the total analog sum of energy 
deposit in the detector. Two signals of more than 500 keV in 200 ^s will be required. 

8.4.7 Summary of the systematic uncertainty cancellations 

A summary of the systematic errors associated with V e event selection cuts is given 
in Table I5~2l We summarize in Table I5~3l the systematic uncertainties that totally 
cancel, or to a large extent, in the Double-CHOOZ experiment. The error on the 
absolute knowledge of the chemical composition of the Gd scintillator disappears. 

1 Energy cut on gamma spectrum from a Gd neutron capture. 



CHOOZ Double-CHOOZ 



selection cut 


rel. error (%) 


rel. error (%) 


Comment 


positron energy* 


0.8 





not used 


positron-geode distance 


0.1 





not used 


neutron capture 


1.0 


0.2 


Cf calibration 


capture energy containment 


0.4 


0.2 


Energy calibration 


neutron-geode distance 


0.1 





not used 


neutron delay 


0.4 


0.1 




positron-neutron distance 


0.3 


0-0.2 


if not used 


neutron multiplicity* 


0.5 





not used 


combined* 


1.5 


0.2-0.3 





* average values 



Table 8.2: Summary of the neutrino selection cut uncertainties. CHOOZ values 
have been taken from CHO03]. 





CHOOZ 


Double-CHOOZ 


Reactor power 


0.7 % 


negligible 


Energy per fission 


0.6 % 


negligible 


V c /fission 


0.2 % 


negligible 


Neutrino cross section 


0.1 % 


negligible 


Number of protons/cm 3 


0.8 % 


0.2 % 


Neutron time capture 


0.4 % 


negligible 


Neutron efficiency 


0.85 % 


0.2 % 


Neutron energy cut 1 


0.4 % 


0.2 % 



Table 8.3: Summary of systematic errors that cancel or are significantly decreased 
in Double-CHOOZ. 



There remains only the measurement error on the volume of target (relative between 
two detectors). The error on the absolute knowledge of the gamma spectrum from 
a Gd neutron capture disappears. However, there will be a calibration error on the 
difference between the 6 MeV energy cut in both detectors. 

8.4.8 Systematic uncertainties outlook 

Table 1^41 summarizes the identified systematic errors that are currently being con- 
sidered for the Double-CHOOZ experiment. 

8.5 Background subtraction error 

The design of the detector will allow a Signal/Background (S/B) ratio of about 
100 to be achieved (compared to 25 at full reactor power in the first experiment 
CHO03 ). The knowledge of the background at a level around 30-50 % will reduce 
the background systematic uncertainties to an acceptable level. In the Double- 
CHOOZ experiment, two background components have been identified, uncorre- 
cted and correlated (see Chapter UJ). Among those backgrounds, one has: 

• The accidental rate, that can be computed from the single event measure- 
ments, for each energy bin. 



After CHOOZ Double-CHOOZ Goal 



Solid angle 


U.Z /0 


to confirm 


1 1 TYl C* 
V Ul LII11L, 


n 9 % 

U.Z /o 


LU L (J 11 11 1 111 


Density 


0.1 % 


0.1 % 


Ratio H/C 


0.1 % 


0.1 % 


Neutron efficiency 


0.2 % 


0.1 % 


Neutron energy 


0.2 % 


0.2 % 


Spatial effects 


neglect? 


to confirm 


Time cut 


0.1 % 


0.1 % 


Dead time(veto) 


0.25 % 


to improve 


Acquisition 


0.1 % 


.1 % 


Distance cut 


0.3 % 


0-0.2 % 



Grand total 0.6 % < 0.6 % (to confirm) 



Table 8.4: The column "After CHOOZ" lists the systematic errors that can be 
achieved without improvement of the CHOOZ published systematic uncertainties 
CHOQlj. In Double-CHOOZ, we estimate the total systematic error on the nor- 
malization between the detectors to be less than 0.5 %. The aim of the work prior 
the final proposal is to confirm this number, and thus increase the safety margin of 
the experiment. 



• The fast neutrons creating recoil protons, and then a neutron capture. This 
background was dominant in the first experiment |CHO03| . The associated 
energy spectrum is relatively flat up to a few tens of MeV. 

• The cosmogenic muon induced events, such as 9 Li and 8 He, that have been 
studied and measured at the NA54 CERN experiment [HagOO in a muon beam 
as well as in the KamLAND experiment KAM02 . Their energy spectrum 
goes well above 8 MeV, and follows a well defined shape. 

The backgrounds that will be measured are: 

• Below 1 MeV (this was not possible in the first experiment, due to the different 
detector design and the higher energy threshold) 

• Above 8 MeV (where there remains only 0.1 % of the neutrino signal). 

• By extrapolating from the various thermal power of the plant (refueling will 
result in two months per year at half power). 

From the measurement of the accidental events energy shape, and from the extrac- 
tion of the cosmogenic events shape, the shape of the fast neutron events can be 
obtained with a precision greater than what is required. 

8.6 Liquid scintillator stability and calibration 

The experiment has some sensitivity to a slight distortion induced by neutrino 
oscillations. A rate only analysis would only provide a sensitivity that is twice the 
quoted value of 0.03 on sin 2 (26>i3). From the simulation, identical energy scales at 
the 1 % level is necessary. The specification of no more than 100 keV scale difference 
at 6 MeV is achieved if this 1 % level is obtained. This relative calibration is 
easier than an absolute linearity, but still very important to consider in the detector 
design. We can, for example, move the same calibration radioactive sources from 
one detector to the other, and directly compare the position of the well defined 
calibration peaks. 



Chapter 9 

Sensitivity and discovery 
potential 

We describe here the details of the simulation of the Doublc-CHOOZ experiment. 
The sensitivity to sin 2 (26*13) is presented in Section l^3l and we present the discovery 
potential of the experiment in Sect ion f9. 41 The statistical analysis (systematic error 
handling) introduced here is based on the work of |Hub02| . 

9.1 The neutrino signal 

In this section we describe the set of parameters used in the simulation. 

9.1.1 Reactor V e spectrum 

The V e spectrum above detection threshold is the result of (3~ decays of 235 U, 
238 U, 239 Pu and 241 Pu fission products. Measurements for 235 U, 239 Pu an d 241 Pu 
and theoretical calculations for 238 U are used to evaluate the v e spectrum |Sch85 
Hah89 . While a nuclear reactor operates, the fission products proportions evolve 
in time; as an approximation in this evaluation, we use a typical averaged fuel 
composition during a reactor cycle corresponding to 55.6 % of 235 U, 32.6 % of 
239 Pu, 7.1 % of 238 U and 4.7 % of 241 Pu. The mean energy release per fission is 
then 203.87 MeV and the energy weighted cross section for V e p — ► ne + amounts to 
(cr) fission = 5.825 • 1CT 43 cm 2 per fission. 

9.1.2 Detector and power station features 

Table l9~T1 contains the principal features of the CHOOZ power station nuclear cores, 
as well as their distances from the near and far detectors. Table l9~2l presents the 





CHOOZ-B-1 


CHOOZ-B-2 


Electrical Power (raw/net GW e ) 


1.516/1.455 


1.516/1.455 


Thermal power (GWth) 


4.2 


4.2 


Global load factor 


80 % 


80 % 


Near detector distance 


100-200 m 


100-200 m 


Far detector distance 


1,000 m 


1,100 m 



Table 9.1: Chooz power station main features |GEA01j . 
characteristics of the detectors used in the simulation. We considered a target 
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scintillator composition of 20 % of PXE and 80 % of dodecane (see Chapter [EJ ■ 
This translates into 8.33 • 10 29 free protons in the 12.7 m 3 inner acrylic vessel. For 
simplicity we assume that the two cores are equivalent to a single core of 8.4 GW t h 
located 150 m away from the near detector and 1,050 m from the far detector. We 
checked that a full simulation with two separated cores at CHOOZ does not change 
the results presented here. The global load factor of the CHOOZ nuclear reactor 





Near Detector 


Far Detector 


Distance 


100 m 


1,050 m 


Target volume 


12.7 m 3 


12.7 m 3 


Target mass 


10.16 tons 


10.16 tons 


Free H 


8.33 10 29 


8.33 10 29 


Detection efficiency 


80 % 


80 % 


Reactor efficiency 


80 % 


80 % 


Dead time 


50 % 


a few % 


Overall efficiency 


32 % 


64 % 


V e events after 3 years 


3,213,000 


58,000 



Table 9.2: Detector parameters used in the simulation. As an example we take here 
the near detector distance at 100 m. Results presented in this chapter don't change 
if this distance is increased to 200 m. 

is taken to be 80 %. We assume that the detection efficiency for both detectors is 
80 % (69.8 % in CHOOZ |CHO03| ). We neglect the dead time for the far detector 
(300 m.w.e. overburden). Since the CHOOZ near site will be shallower, between 
60 to 80 m.w.e, we apply a dead time of 50 % to be conservative (a 500 /isec cut 
to each muon crossing the detector leads to a dead time around 30 % at 60 m.w.e). 
The overall efficiencies used in the simulation for the near and far detectors are thus 
respectively 32 % and 64 %. 

9.1.3 Expected number of events 

Neglecting the correction terms of order a = (^ £™S° l 1 ~ (2 • 10 -2 ) 2 , we used the 
following V c survival probability: 

9/ , , / ATOnoLlml \ 
= 1 ~ Sin (2013) Sin ( L27 l#eV ) ' (9 ' 1} 

The expected number of antineutrino events in the near (Nf) and far detector 
(Nf), in the energy bin [Ei,E i+ {\, is 

N A =T A / S(E U , E' v )a(E u )^(E u , L a )P v ,^ Vq (E u ,L a ) &E v &K, (9-2) 

JEi JO 

where A — N,F. The cross section a is given in equation 12. Ill and the V e flux is 
computed according to Figures 12.11 and 12.21 . The normalization factor T includes 
the global load factor G (fraction of running time of the reactors over a year), the 
reactor thermal power P, the detector efficiency e A , the dead time fraction D A , 
the target volume V and the exposure time T: 

F A = GxPxVxTx(l- D A ) x e A (9.3) 

The energy resolution effect is taken into account as follows: 

S{E,E') = AT ( E - E\?-p) , (9.4) 



where Af is a Gaussian distribution. In practice, we have used an energy bin size at 
least four times larger than the energy resolution effect and thus we neglected it in 
first approximation for this analysis. We checked this approximation by comparing 
our r esults with the work of |Hub02llHiiE04| . 



9.2 Systematic errors handling 
9.2.1 x 2 analysis 

In this section we describe the x 2 -analysis of the near-far detector set and how we 
implemented the systematic errors previously discussed. We write Of the computed 
number of events observed in i th energy bin in near (A — N) and far (A = F) 
detectors. The theoretical predictions for the detector A in the z th bin is 

Scores 



T t A = (1 + a + b A + a) £ (1 + h)N A + g A M{ 

3 = 1 



(9.5) 



where a, b , c,-, fj, g A will be the fitted parameters. M A is the first order correction 
term to take into account the energy scale uncertainty obtained by replacing E w - ls 

by ;i • g A )E viB : 



N, 



Mi 



v . dNA(g A ) 



dg J 



(9.6) 



where Nfj is the computed number of events in i bin in detector A coming from 
the j th reactor core: 



N A - 



E i+1 f+oo 



S(E u ,El)a(E u )^(E^ L A )P Ve ^ [E v , L A ) AE v dE' v 



(9.7) 



which depends on the oscillation parameters through the survival probability. The 
observed number of events Of has been chosen to be the computed event number for 
given "true values" of the oscillation parameters: Of = Y^j^T" Njfjisvn 2 (26*13), Am 2 ). 
We used a x 2 function including the full spectral information from both detectors: 
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(9.8) 



For each point in the oscillation parameters space, the \ 2 function has to be 
minimized with respect to the parameters a, b N , b F , a, g N , g F , d, e N , e F , fj 
modeling the systematic errors. The parameter a refers to the error on the overall 
normalization of the number of events common to both detectors. Parameters b N 
and b F relate to the uncorrelated normalization uncertainties of the two detectors. 
The energy scale uncertainty is taken into account through parameters g N and 



g F in the expression of in equation 19.51 We assumed here a flat background 
distribution: 



The numerical minimization has been performed with the MINUIT package |.1 am| . 
We now discuss all the relevant terms of Eauation l9.9l in turn. 

9.2.2 Absolute normalization error: cr a b s 

We include a common overall normalization error for the event rate of the near 
and far detectors. This error accounts for the uncertainty on the v e flux of the 
reactor, the detection cross section, or any bias that could affect both detectors in 
the same way 1 . This error is of the order of a few percent; one has for instance 1.4 % 
in |Dec94| . 2 % in [KAM02 . The overall normalization error has almost no impact 
on the sensitivity to an oscillation effect in the Double-CHOOZ experiment since 
two detectors will be used (see Figures 19.71 and l9.8J) . Nevertheless, we included an 
absolute normalization error <7 a b s = 2 % in the simulation. 

9.2.3 Relative normalization error: (x rel 

We take into account an uncorrelated normalization uncertainty between the near 
and far detectors. This is the dominant experimental error for the Double-CHOOZ 
experiment. There are contributions from uncertainties on the detector design (fidu- 
cial volume, stability of the liquid scintillator, relative dead time measurement) and 
the uncertainties related to the v e event selection cuts (relative detector efficiency) . 
According to the results presented in Chapter [5] we take the relative normalization 
error a Ic i — 0.6 % as our default value. 

9.2.4 Spectral shape error: cr shp 

To take into account the v spectrum shape uncertainty, we introduce an error a s hp 
on the theoretical prediction for each energy bin which we take to be fully uncor- 
related between different energy bins. Since this error is induced by the physical 
uncertainty on the fission product beta decay spectra, it is fully correlated between 
the corresponding bins in the near and far detector. In the simulation we use the 
shape error value a s \ lp — 2 %, as measured in |BUG96 . 

9.2.5 Energy scale error: ct sc i 

We take into account the energy scale calibration uncertainty by introducing a 
parameter g A for each detector (A = N,F), and replacing the observed energy 
E bs by (1 + ,g j4 )-E bs- We assume that the energy calibration is known with an 
error of ct sc i ~ 0.5 %. We found that, as long as no detailed background simulation 
is performed on the data, this error can be neglected in first approximation for 
the sensitivity computations. This is understandable since the Double-CHOOZ 
experiment is mostly sensitive to the number of events integrated over the whole 
positron spectrum. Nevertheless, a careful study of this error is going on to better 
understand its influence on the discovery potential of Double-CHOOZ. 

1 For instance, a bias in the volume measurement affecting the two detectors is equivalent to an 
uncertainty in the reactor V e flux. 




(9.9) 



9.2.6 Individual core power fluctuation error: a C R 



Since the Double-CHOOZ power station has two nuclear cores, we introduced an 
independent error of a c a = 0.5 % mimicking a thermal power fluctuation of each 
nuclear core. Indeed, depending on the exact location of the near detector site, 
the near and far detectors will not receive the same V c contribution from both 
cores. In that case, an independent fluctuation in the two cores could lead to a 
relative systematic error between the detectors. However, we found this error to be 
negligible and we do not consider it further. 



9.2.7 Background subtraction error 

We considered two different ways to introduce an error on the background subtrac- 
tion procedure. 

9.2.7.1 Reactor V e shape background: Ob2b 

This is modeled as an uncorrelated error <Jb2b in the background subtraction step. 
This error is bin-to-bin uncorrelated, uncorrelated between the near and far detec- 
tors, and proportional to the bin content (i.e. the background has the same shape 
as the positron spectrum). Typically we used values ranging from Ob2b = 0.5 % to 
Cb2b = 1-5 %. 



9.2.7.2 Flat background: Obkg 

This background is closer in shape to the background of fast neutrons created in the 
rocks close to the detector. It was dominating in CHOOZ CHO03 , and is expected 
to play an important role in Double-CHOOZ as well (see Chapter EJ). We assume 
that it amounts typically for Rjy = 1 % & Rp = 1 % of the total V e signal. To be 
conservative we consider an error on those rates of cr^ kg = 100 % a£ kg = 100 %, in 
the near and far detectors. A careful study of the impact of the background on the 
sensitivity and on the discovery potential as well is going on. 



9.3 Sensitivity in the case of no oscillations 

We present our results for the current best fit value of the atmospheric mass splitting 
Am 2 ^ = 2.01q 7 ' 10~ 3 °V 2 _SK02a as our default value. Nevertheless, we also used 
the recent analysis of the Super-Kamiokande data leading to Am| 3 = 2. 4^0^ ■ 10 _3 eV 2 
SK04a , for completeness. We also assume that a forthcoming accelerator exper- 
iment will provide a precise measurement of Ato 2 3 , with an error better than 
20 %: a Am 2 = 0.2 • Am 2 prior to the Double-CHOOZ result jMINOlal IMINOlbj . 
Figure 19.11 displays the expected sensitivity of Double-CHOOZ in the case of no- 
oscillations, as a function of time. In this case we have a sensitivity of sin 2 (2#i3) < 
0.045 (90 % C.L.) after one year of data taking, and sin 2 (2#i 3 ) < 0.03 after three 
years. The sensitivity dependence with respect to the atmospheric mass splitting 
adm2 value is shown in Figure 19.31 Figure 19.41 displays the effect of a re i on the 
sensitivity of Double-CHOOZ in the (sin 2 (26>i3), Am 2 tm ) plane. The relative nor- 
malization influence on the sin 2 (2^13) limit as a function of the exposure time is 
shown in Figure 1^31 









CM CM ( 

E E 

< < 


= 2.0 10 

= 1.5 10 


1 eV 2 J 

:ev: ■ 


\ 












= 3.0 1 C 


" J eV • 



















































































































































1 23456789 10 

Exposure time in years 



Figure 9.1: Evolution of sin 2 (2f5i3) sensitivity with the exposure time. The three 
curves shown here are for different values of Am 2 as shown in the legend. 
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Figure 9.2: Evolution of sin 2 (2f5 13 ) sensitivity with the exposure time. The three 
curves shown here are for different values of Am 2 as shown in the legend. These val- 
ues have been chosen from the second analysis (L/E) of the same Super-Kamiokande 
data |SK04aj . 



9.3.1 Comparison of Double-CHOOZ and the T2K sensitiv- 
ities 

We compute both the Double-CHOOZ and the T2K sensitivities, in the sin 2 (2f3 13 )-c) 
plane, for three dates: January 2009, January 2011, and January 2015. We assume 
that the Double-CHOOZ experiment will start to take data with two detectors on 
January 2008, while the T2K experiment will start exactly two years later, on Jan- 
uary 2010, with the nominal beam intensity (since the T2K neutrino line is expected 
to be completed within the year 2009, we assume that the accelerator commission- 
ing will be finished by the end of 2009 |SK04bl [ST^04c| ). For the computation of 
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Figure 9.3: Double-CHOOZ sensitivity limit at 90 % C.L. (for 1 d.oi). 
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Figure 9.4: Influence of the relative normalization uncertainty on the sin 2 (26*13) 
limit in the (sin 2 (26*13), Am 2 ) plane in the case of no oscillations (for three years of 
operation). 



the Double-CHOOZ sensitivity we assume here a relative normalization error of 
0.6 % for both detectors. The correlated backgrounds considered here amount to 
1.5 % of the signal for both the near and far detectors. Several background com- 
ponents of known shape have been included (proton recoil, accidental, spallation, 
see Chapter 0). An additional uncorrelated background component of 0.5 % is also 
considered here. All backgrounds are supposed to be known with a 50 % error. De- 
tails of the analysis procedure are given in |Hub021 lHub04| . For the simulation of 
the T2K experiment, the experimental parameters are taken from |T2K02I IT2K03 . 
We used nominal 1 year and 5 year running times for T2K, and 1, 3, and 7 years 
for the reactor setup (with 20,000 events/year). We compute the two-dimensional 
allowed fit regions (i.e., the parameters on the axes are the fitted parameters, in the 
sin 2 (26'i3)-(5 plane) for three dates: January 2009, January 2011, and January 2015. 
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Figure 9.5: Influence of the relative normalization uncertainty on the sin (26*13) 
limit as a function of the exposure time (in years) in the case of no oscillations. 



The curves for T2K include all correlations and degeneracies and are obtained as 
projections of the fit manifolds onto the sin 2 (26* 13 )-r5 plane Hub02i lHub04| . 
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Figure 9. 6: Limi t at 90 % C.L. in the sin 2 (26 13 )-6 plane for Double-CHOOZ 
and T2K |Hub02l lHub04j . The following oscillation parameters have been used: 
Am§! = 2 • 10" 3 eV 2 , Am^ = 7 • 10" 5 eV 2 , sin 2 (2(9 23 ) = 1.0, sin 2 (26>i 2 ) = 0.8, 



and sin (26*13) = 0. We have considered 1 d.o.f for the analysis of the Double- 
CHOOZ experiment, but 2 d.o.f. for the analysis of T2K that is sensitive to both 
sin 2 (26*13) and 5 simultaneously. 90 % C.L. intervals are shown with solid lines, 
and 3cr intervals are displayed with dashed lines. The thick curves describe the 
Double-CHOOZ setup, and the thin curves the T2K experiment, with black curves 
for best-fit solution, and gray curves for the sgn(Am§i)-degeneracy 



9.4 Discovery potential 

9.4.1 Impact of the errors on the discovery potential 

The 3cr discovery potential of Double-CHOOZ is displayed on Figures l£T7l and fD~51 

for respectively Am§i = 2.0 and2.4 • 10 -3 eV 2 . In the first case, a non- vanishing 



value of sin 2 (26*13) = 0.05 could be detected at 3 a after three years of data taking. 
For the second case, this value becomes sin 2 (26*i 3 ) = 0.04. 

9.4.2 Comparison of Double-CHOOZ and the T2K discovery 
potential 

The computation is done as presented in Section ®. 3. II for both the Double-CHOOZ 
and the T2K experiments taken at three dates: January 2009, January 2011, and 
January 2015. To investigate the discovery potential of both experiments, we used 
three benchmark values sin 2 (26*13) = 0.14, 0.08, 0.04. Results are presented respec- 
tively in the sin 2 (20i 3 )-<5 plan in Figures EHUl and ETTT1 
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Figure 9.7: Statistical and systematic errors contributions to sin 2 (2#i 3 ) measure- 
ment. We assumed here SK-I analysis best fit value Am^ = 2.0 10~ 3 eV 2 , 3 years 
of data taking for Double-CHOOZ with 64 % (expecting around 58,000 events in 
the case of no oscillations) of efficiency in the far detector and 32 % in the near one. 
We also set the systematic errors to the standard ones: the absolute normalization 
to 2 %, the relative to 0.6 %, the shape uncertainty to 2 % and the background to 
1 %. The different error intervals are plotted at with a 3 a confidence level. We see 
here that the discovery potential limit of Double-CHOOZ to detect a non- vanishing 
value of sin 2 (2^13) is around 0.05. We also see here that struggling harder than 
the level of 0.6 % on the relative normalization could lower this discovery potential 
limit. 
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Figure 9.8: Same as Figure l9~7l but for Ar 
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Figure 9.9: Measurement of sin 2 (20i 3 ) and 5 with Double-CHOOZ and T2K Hub02 
IHub04| . The following oscillation parameters have been used: Am^j = 2 • 10 -3 eV 2 , 
Amii = 7- 10" 5 eV 2 , sin 2 (2(9 23 ) = 1.0, sin 2 (26>i 2 ) = 0.8. The 13 mixing angle was 
generated as sin 2 (2#i3) = 0.14 and the CPS phase has been fixed at S = tt/2. We 
considered 1 d.o.f. for the analysis of the Double-CHOOZ experiment, but 2 d.o.f. 
for the analysis of T2K that is sensitive to both sin 2 (2^13) & S simultaneously. 90 % 
C.L. interval are shown with solid lines, and 3<r intervals are displayed with dashed 
lines. The thick curves describe the Double-CHOOZ setup, and the thin curves the 
T2K experiment, with black curves for best-fit solution, and gray curves for the 
sgn(Am| 1 )-degeneracy. The minimum \ 2 is drawn at marked points. 
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Figure 9.10: Same as Figure 
sin 2 (26»i 3 ) = 0.08. 
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Figure 9.11: Same as Figure l9~iH but for the #13 mixing angle was generated at 
sin 2 (26»i 3 ) = 0.04. 



Appendix A 

z7 e and safeguards 
applications 

The International Atomic Energy Agency (IAEA) is the United Nations agency in 
charge of the development of peaceful use of atomic energy IAEA . In particular 
IAEA is the verification authority of the Treaty on the Non-Proliferation of Nuclear 
Weapons (NPT). To do that job inspections of civil nuclear installations and related 
facilities under safeguards agreements are made in more than 140 states. IAEA use 
many different tools for these verifications, like neutron monitors, gamma spec- 
troscopy, but also bookkeeping of the isotopic composition at the fuel element level 
before and after their use in the nuclear power station. In particular it verifies that 
weapon-origin and other fissile materials that Russia and USA have released from 
their defense programs are used for civil applications. 

The existence of a v c signal sensitive to the power and isotopic composition of a 
reactor core could provide a mean to address certain safeguards applications. Thus 
the IAEA very recently asked member states to make a feasibility study to deter- 
mine whether antineutrino detection methods might provide practical safeguards 
tools for selected applications. If this method proves to be useful, IAEA has the 
power to decide that any new nuclear power plant to be built has to include an V e 
monitor. 

The high penetration power of antineutrinos and the detection capability might 
provide a mean to make "remote" and non- intrusive measurements of plutonium 
content in reactors and in large inventories of spent fuel. The antineutrino flux and 
energy spectrum depend upon the thermal power and the fissile isotopic composi- 
tion of the reactor fuel. Because the antineutrino signal from the reactor decreases 
as the square of the distance from the reactor to the detector the " remote" measure- 
ment is really only practical at distances of a few tens of meters if one is constrained 
to "small" detectors of the order of few cubic meters in size. Based on predicted 
and observed spectra, the number of v e per fission from 239 Pu is known to be less 
than the number from 235 U. This variation has been directly measured in reactor 
antineutrino experiments. This may offer a mean to monitor changes in the rela- 
tive amounts of 235 U and 239 Pu in the core and in freshly discharged spent fuel. If 
made in conjunction with accurate independent measurements of the thermal power 
(including the ambient reactor temperature and the flow rate of cooling water), an- 
tineutrino measurements might provide an estimate of the isotopic composition of 
the core, in particular plutonium inventories. The shape of the antineutrino spec- 
trum can provide additional information about core fissile isotopic composition. 
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In order to determine the feasibility of antineutrino detection for safeguards ap- 
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241p u 


v c /fission 
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9.0 
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9.3 



Table A.l: Number of v c emitted per fission and end points of U and Pu fissile 
isotopes. 

plications, a series of scenarios involving antineutrino detectors should be defined, 
both for reactors and for spent fuel inventories. The effectiveness, sensitivity, and 
possible vulnerabilities of antineutrino detection should be examined for these sce- 
narios. For the IAEA, the proposed feasibility study should seek to establish or 
refute the utility of antineutrino detection methods as a new safeguards tool, and 
serve as a guide for future efforts. Additional lab tests and theoretical calculations 
should also be performed to more precisely estimate the underlying j3 spectra of 
plutonium and uranium fission products, especially at low energies, corresponding 
to the most energetic antineutrinos. 

The appropriate starting point for this scenario is a representative PWR. For 
this reactor type, simulations of the evolution of the antineutrino flux and spectrum 
over time should be provided, and the required precision of the antineutrino detec- 
tor and independent power measurements should be estimated. In that respect the 
measurement performed by the Double-CHOOZ experiment with its near detector, 
as it is explained in the proposal, will constitute the most precise determination of 
the antineutrinos emitted by a PWR. In particular, the follow-up of the spectrum 
and rate after refueling with fresh 235 U, would allow a precision study of the corre- 
lation between plutonium content and the measured spectrum. If it is possible in 
addition to have a detailed follow-up of the evolution of the fuel burn-up, by the 
use of fission chambers, the data gathered by these experiments will constitute an 
excellent experimental basis for the above feasibility studies of potential monitor- 
ing and for bench-marking fuel management codes. This measurement will help to 
meet another important point of the IAEA concern, linked to the verification of 
provisions of the US-Russian Plutonium Management and Disposition Agreement 
(PMDA). This agreement concerns MOX fuel made using weapon origin plutonium. 

Verifying core burn up while the reactors are operating would provide a mean to 
determine whether or not the disposition criteria have been met. From the present 
knowledge of the antineutrino spectrum emitted by the fission products, we see that 
the most energetic part offers the best possibility to disentangle fission from 235 U 
and 239 Pu. Unfortunately the present uncertainty in that region of energy is rather 
large, due to the difficulties of measuring the corresponding low energy /3 _ . 

Thus, in relation to this feasibility studies, new measurements of the /3 spec- 
trum for the various fissile elements are mandatory. A group of nuclear physicists 
has developed tools, in the frame of MinilNCA collaboration |Inca| . which can be 
modified to perform these measurements at ILL. Needless to say that a more pre- 
cise knowledge of the antineutrinos emitted in the reactor core would also benefit 
the physics measurements of #13. The overall IAEA feasibility studies are larger 
than the topics briefly described above. It is also of interest to study other present 
reactor types, like BWRs, FBRs, and possibly CANDU reactors. Future reactors 
(e.g., PBMRs, Gen IV reactors, accelerator-driven sub-critical assemblies for trans- 
mutation), especially reactors using carbide, nitride, metal or molten salt fuels must 
also be considered. IAEA seeks also to the possibility of monitoring large spent- 



fuel elements. For this application, the likelihood is that antineutrino detectors 
could only make measurements on large quantities of emitters, e.g., several cores 
worth of spent fuel. In the time of the experiment the discharge of parts of the 
core will happen and the Double-CHOOZ experiment will quantify the sensitivity 
of such monitoring. More generally the techniques developed for the detection of 
antineutrinos could be applied for the monitoring of nuclear activities at the level 
of a country. For example a KamLAND type detector |KAM02| deeply submerged 
off the coast of the country, would offer the sensitivity to detect a new underground 
reactor located at several hundreds of kilometers. In that respect, the progress in 
term of detecting medias (Gd doped liquid scintillators) would be greatly helpful. 
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Nuclear reactor j3 spectra 



New measurements of the (3 spectrum for various fissile elements present in a nuclear 
reactor will be very important for the Double-CHOOZ experiment to understand 
the physics at the near detector. Of course, it is less important for the oscillation 
analysis, since the absolute normalization error on the V c flux is absorbed if two 
detectors are used simultaneously at different baselines. These new integral mea- 
surements deal with a complete characterization of the (3 spectrum produce in the 
fuel element by taking into account the evolution of the fuel. This information is im- 
portant to characterize the antineutrino spectra at the Double-CHOOZ experiment 
but is also unavoidable for the feasibility studies of using antineutrino detection 
methods as a new safeguards tool. 

In the frame of the Mini-INCA project |lnca| . the group has developed a set 
of experimental tools to perform quasi online a- and 7-spectroscopy analyzes on 
irradiated isotopes and to monitor online the neutron flux in the high flux reactor of 
the ILL reactor. It has also developed competences on the Monte-Carlo simulations 
of complex systems and in particular nuclear reactors. These competences will be 
used to provide to the community a set of integral f3 energy spectra relevant for 
the Double-CHOOZ experiment and for safeguards studies and to understand and 
monitor all the fluctuations in the antineutrino spectra originated from the reactor 
source. 



B.l New (3 energy spectra measurements at ILL 

The a and 7 spectroscopy station, connected to an irradiation channel of the ILL 
reactor, offer the possibility to perform irradiations in a quasi thermal neutron flux 
up to 20 times the nominal value in a PWR. This irradiation can be followed by 
measurements and repeated as many time as needed. It offers then the unique 
possibility to characterize the evolution of the beta spectrum as a function of the 
irradiation time and the irradiation cooling. The expected modification of the (3 
spectrum as a function of the irradiation time is connected to the transmutation 
induced by neutron capture of the fissile and fission fragment elements. It is thus 
related to the natural evolution of the spent-fuel in the reactor. The modification 
of the (3 spectrum as a function of the cooling time is connected to the decay chain 
of the fission products and is then a mean to select the emitted fragments by their 
time of live. This information is important because long-lived fission fragments 
accumulate in the core and after few days mainly contribute to the low energy 
part of the antineutrino-spectra. We propose to modify the spectroscopy station by 
adding a large dynamic (3~ spectrometer and to measure the (3 spectra for 235 U, 
239 Pu, 241 Pu and 243 Cm for different irradiation and cooling times. Due to the 
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mechanical transfer of the sample from the irradiation spot to the measurement 
station an irreducible delay time of 30 mn is imposed leading to the loss of short- 
live fragments. To characterize the (3 prompt emissions online measurements will 
be done on a neutron guide where cold neutrons are available. 

B.2 Reactivity monitoring 

Micro-fission chambers developed for high neutron fluxes are used in core in the ILL 
reactor. They provide very precise neutron flux measurements and allow to monitor 
in line the reactivity fluctuations of the core. Due to their small dimensions (4 mm 
in diameter and 4 cm in length) and the low fissile deposit, they should allow 
to measure very precisely the gravity center of the core, with a negligible flux 
perturbation, if placed out core of the Chooz reactor. 

B.3 Double-CHOOZ reactor core simulation and 
follow-up 

By the mean of Monte-Carlo and deterministic codes developed for neutron flux 
calculation and evolution at ILL and for various type of transmutation scenario, we 
propose to model the complete history of Chooz reactor core to study the sensitivity 
of the neutrino spectrum to the isotopic composition and fuel burn up. 



Appendix C 



Some numbers from the 
CHOOZ experiment 

The CHOOZ experiment |CHQ98I ICHQ99I ICHO00I ICHO03j was located close to 
the CHOOZ nuclear power plant, in the North of France, 10 km from the Belgian 
border. The power plant consists of two twin pressurized water reactors (PWR), 
the first of a series of the newly developed N4 PWR generation in France CEA01 . 
The thermal power of each reactor is 4.25 GW C (1.3 GW e ). These reactors started 
respectively in May and August 1997, just after the start of the data taking of the 
CHOOZ detector (April 1997). This opportunity allowed a measurement of the 
reactor-off background, and a separation of individual reactors contributions. 

The detector was located in an underground laboratory about 1 km from the 
neutrino source. The 300 m.w.e. rock overburden reduced the external cosmic 
ray muon flux, by a factor of about 300, to a value of 0.4 m~ 2 s _1 . This was the 
main criterion to select this site. Indeed, the previous experiment at the Bugey 
reactor power plant |BUG96| showed the requirement of reducing by two orders of 
magnitude the flux of fast neutrons produced by muon-induced nuclear spallations in 
the material surrounding the detector. The neutron flux was measured at energies 
greater than 8 MeV and found to be about 1/day, in good agreement with the 
prediction. 

The detector envelope consisted of a cylindrical steel vessel, 5.5 m diameter and 
5.5 m height. The vessel was placed in a pit (7 m diameter and 7 m deep), and 
was surrounded by 75 cm of low activity sand. It was composed of three concentric 
regions, from inside to outside: 

• a central 5 tons target in a transparent Plexiglas container filled with a 0.09 % 
Gd-loaded scintillator 

• an intermediate 70 cm thick region, filled with non-loaded scintillator and used 
to protect the target from PMT radioactivity and to contain the gammas from 
neutron capture on Gd. These 2 regions were viewed by 192 PMTs 

• an outer veto, filled with the same scintillator. 

The scintillator showed a degradation of the transparency over time, which re- 
sulted in a decrease of the light yield (live time around 250 days). The event position 
was reconstructed by fitting the charge balance, with a typical precision of 10 cm 
for the positron and 20 cm for the neutron. Source and laser calibrations found that 
due to the small size of the detector the time reconstruction was less precise than 
expected. The reconstruction became more difficult when the event was located 
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near the PMTs, due to the 1/r 2 divergence of the light collected (see Figure 31 of 
The final event selection used the following cuts: 

• positron energy smaller than 8 MeV (only 0.05 % of the positrons have a 
higher energy) 

• neutron energy between 6 and 12 MeV 

• distance from the PMT support structure larger than 30 cm for both positron 
and neutron 

• distance between positron and neutron smaller than 100 cm 

• low particles multiplicity: when a third particle is detected in the time window 
between the positron and neutron candidates, a complicated cut must be 
applied (see 8.7 of |CH()03p . 

The neutron capture on Gd is identified by a 6 MeV cut on the total energy 
emitted. This cut induce a systematic error of 0.4 %, due to the poor knowledge of 
the emission spectrum of the gammas released after the neutron capture. 

The scintillating buffer around the target was important enough to reduce the 
gammas escape. This cut was calibrated with a neutron source. The 3 cuts on the 
distances were rather difficult to calibrate, due to the the reconstruction problems 
described above. This created a tail of badly reconstructed events, which was very 
difficult to simulate (0.4 % systematic error on the positron-neutron distance cut). 
The positron threshold was carefully calibrated, as shown in Figure 39 of [UHO03 . 
The value of the threshold depends upon the position of the event, due to the 
variation of solid angle and to the shadow of some mechanical pieces such as the 
neck of the detector (0.8 % systematic error). The time cut relied on Monte-Carlo 
simulation. The corresponding systematic error was estimated to be 0.4 %. The 
final result was given as the ratio of the number of measured events versus the 
number of expected events, averaged on the energy spectrum. It was found to be: 

R = 1.0 ± 2.8 % (stat) ± 2.7 % (sys). 

Two components were identified in the background: 

• Correlated events: which had a flat distribution for energies greater than 
8 MeV, and were due to the recoil protons from fast spallation neutrons. It 
was extrapolated to 1 event/day. 

• Accidental events: which were obtained from the measure of the singles rates. 

The total noise was measured during the reactor-off, and by extrapolating the signal 
versus power straight line (see Figure 49 of |CHO03p . It is in good agreement with 
the sum of the correlated and accidental components. These numbers have to be 
compared to a signal of 26 events/day at full reactor power. The systematic error 
was due mainly to the reactor uncertainties (2 %), the detector efficiency (1.5 %), 
and to the normalisation of the detector dominated by the error on the proton 
number from the H/C ratio in the liquid (0.8 %). The resulting exclusion plot is 
shown in Figure 58 of |CH(J03| . The corresponding limit on sin 2 (2^13) is 0.14 for 
Am 2 = 2.6 10~ 3 eV 2 , and 0.2 for Am 2 = 2.0 10~ 3 cV 2 . This limit disappears for 
Am 2 < 0.8 10~ 3 cV 2 , due to the ~1 km distance between the cores and the CHOOZ 
detector. 
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